Tau Functions and Matrix Integrals 

A. Yu. Orlov* 



(N 

o 
o 

(N 
> 

Q 

o 

(N 



Abstract 

We consider solvable matrix models. We generalize Harish-Chandra-Itzykson-Zuber 
and certain other integrals ( Gross- Witten integral and integrals over complex matrices) 
using the notion of tau function of matrix argument. In this case one can reduce matrix 
integral to the integral over eigenvalues, which in turn is certain tau function. The result- 
ing tau functions may be analyzed either by the method of orthogonal polynomials or by 
Schur functions expansion method. 
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1 Introduction 

This paper in its main part is based on the previous paper " Soliton theory. Symmetric Functions 
and Matrix Integrals" |1[ and the papers 0,0. In the first quoted paper we considered certain 
class of tau-functions which we call tau-functions of hypergeometric type. The key point of 
the paper was to identify scalar products on the space of symmetric functions with vacuum 
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expectations of fermionic fields. We considered some known matrix models, hypergeometric 
functions of matrix argument, Schur function expansion as the perturbation series for matrix 
integrals. 

We enlarge the class of matrix integrals which turned out to be tau functions. (This presen- 
tation in particular means that one has a determinant representation, Schur function expansion 
for perturbation series and the possibility to apply the orthogonal polynomials method for non- 
perturbative analysis.) 

Here we include a version of the previous paper into the text to make it more self 
consistent. 

Symmetric functions. Symmetric function is a function of variables x = Xi, . . . , which 
is invariant under the action of permutation group Sn- Symmetric polynomials of n variables 
form a ring denoted by A„. The number n is usually irrelevant and may be infinite. The ring 
of symmetric functions in infinitely many variables is denoted by A. 

Power sums are defined as 

Pm = E^r, m = l,2,... (1.0.1) 

i=l 

Functions PmiPm2 " " " Pm^, k = 1,2, . . . form a basis in A. 

There are different well investigated polynomials of many variables such as Schur functions, 
projective Schur functions, complete symmetric functions Each set form a basis in A. 

There exists a scalar product in A such that Schur functions are orthonormal functions: 

< sx,s^>= 6x,^ (1.0.2) 

The notion of scalar product is very important in the theory of symmetric functions. The 
choice of scalar product gives different orthogonal and bi-orthogonal systems of polynomials. 

Soliton theory. KP hierarchy of integrable equations, which is the most popular example, 
consists of semi-infinite set of nonlinear partial differential equations 

dt^u = K^[u], m = l,2,... (1.0.3) 

which are commuting flows: 

[dt,,d,Ju = (1.0.4) 
The first nontrivial one is Kadomtsev-Petviashvili equation 

dt,u = \dlu + Id^.'dlu + ^3,^ (1.0.5) 



which originally served in plasma physics [|T7|. This equation, which was integrated in 



now plays a very important role both, in physics (see JTSf; see review in |jT^ for modern 



applications) and in mathematics. Each members of the KP hierarchy is compatible with each 
other one. 

Another very important equation is the equation of two-dimensional Toda lattice (TL) 

dtAl^n = e'^"-'-'^" - e'^"-'^"+i (1.0.6) 

This equation gives rise to TL hierarchy which contains derivatives with respect to higher times 
t\,t2, ■ ■ ■ and t*, • • •• 

The key point of soliton theory is the notion of tau function, introduced by Sato (for KP 
tau- function see [0]). Tau function is a sort of potential which gives rise both to TL hierarchy 
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and KP hierarchy. It depends on two semi- infinite sets of higher times ti, ^2, • • • and t^, • • 
and discrete variable n: t = r(n, t,t*). More exphcitly we have [20[],p7|: 



u 



2dl logr(n,t,t* 



/>n(t,t* 



log 



r(n + l,t,t*: 



;i.o.7) 



r(n,t,t*) 

In soliton theory so-called Hirota-Miwa variables x, y, which are related to higher times as 

mt^ = J2xT, mt*^ = J2yr, (1.0.8) 



are well-known ||2J]. Tau function is a symmetric function in Hirota-Miwa variables, higher 
times tm,tl^ are basically power sums, see ( |1.0.1| ), see |^ as an example. 

It is known fact that typical tau function may be presented in the form of double series in 
Schur functions [EHl E9l BOH: 



rfn, t, t* 



;i.o.9) 



where coefficients Kxfj, are independent of higher times and solve very special bilinear equations, 
equivalent to Hirota bilinear equations. 

In the previous paper we use these equations to construct hypergeometric functions which 
depend on many variables, these variables are KP and Toda lattice higher times. Here we 
shall use the general approach to integrable hierarchies of Kyoto school M, see also 



Especially a set of papers about Toda lattice [^, ^ |3T|, ^ Q is important for us. Let 

us notice that there is an interesting set of papers of M.Adler and P. van Moerbeke, where links 
between solitons and random matrix theory were worked out from a different point of view (see 
|l3l,|T^ and other papers). 

Random matrices. In many problems in physics and mathematics one use probability 
distribution on different sets of matrices. We are interested mainly in two matrix models. The 
integration measure dMi 2 is different for different ensembles 



where 



;i.o.io) 

1.0.11) 



Fi(Mi) = Tr ^ M^tm. V2{M2) = Tr ^ M^^t*^ 

m=l m=l 

while U may be different for different ensembles. 

The crucial point when evaluating the integral ( |1.0.10|) is the possibility to reduce it to the 
integral over eigenvalues of matrices Mi, M2. Then ( |1.0.10| ) takes the form 



Z = C 



i=l 



1.0.12) 



In the present paper we are interested in the perturbation series of different matrix models 
which are series in parameters tm^t*^- 



2 Symmetric functions 

2.1 Schur functions and scalar product [j^ 

Partitions. Polynomial functions of many variables are parameterized by partitions. A parti- 
tion is any (finite or infinite) sequence of non-negative integers in decreasing order: 

A = (Ai,A2,...,A^,...), Ai > A2 > ... > A, > ... (2.1.1) 
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Non-zero Aj in ( ^.I.ID are called the parts of A. The number of parts is the length of A, denoted 
by /(A). The sum of the parts is the weight of A, denoted by |A|. If n = |A| we say that A is the 
partition of n. The partition of zero is denoted by 0. 

The diagram of a partition (or Young diagram) may be defined as the set of points (or 
nodes) (i,j) e Z'^ such that 1 < j < Aj. Thus Young diagram is viewed as a subset of entries 
in a matrix with 1{X) lines and Ai rows. We shall denote the diagram of A by the same symbol 
A. 

For example 



U 

is the diagram of (3, 3, 1). The weight of this partition is 7, the length is equal to 3. 

The conjugate of a partition A is the partition A' whose diagram is the transpose of the 
diagram A, i.e. the diagram obtained by reflection in the main diagonal. The partition (3, 2, 2) 
is conjugated to (3, 3, 1), the corresponding diagram is 



There are different notations for partitions. For instance sometimes it is convenient to 
indicate the number of times each integer occurs as a part: 

A = (l'"^2'^^..r"^'■...) (2.1.2) 

means that exactly of the parts of A are equal to i. For instance partition (3, 3, 1) is denoted 
by (1^3^). 

Another notation is due to Frobenius. Suppose that the main diagonal of the diagram of A 
consists of r nodes (i, i) (1 < i < r). Let = Xi—i be the number of nodes in the ith row of A 
to the right of {i, i), for 1 < i < r, and let (3i = X[ — i be the number of nodes in the ith column 
of A below (i, i), for 1 < i < r. We have ai > 0:2 > ■ ■ ■ > > and /3i > > ■ ■ ■ > A' > 0- 
Then we denote the partition A by 

A = (ai,...,a,|/3i,...,/3,) = (a|/3) (2.1.3) 

One may say that Frobenius notation corresponds to hook decomposition of diagram of A, where 
the biggest hook is (ai|/5i), then (a2|/32), and so on up to the smallest one which is {ar\Pr)- 
The corners of the hooks are situated on the main diagonal of the diagram. For instance the 
partition (3, 3, 1) consists of two hooks (2, 2) and (0, 1): 



I I and I I I 

U 

In Frobenius notation this is (2, 0|2, 1). 

Schur functions, complete symmetric functions and power sums. We consider 
polynomial symmetric functions of variables Xi,i = 1,2, ... ,n , where the number n is irrelevant 
and may be infinity. The collection of variables Xi, X2, . . . we denote by x^. When N is irrelevant 
we denote it by x. 
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The symmetric polynomials with rational integer coefficients in n variables form a ring which 
is denoted by A„. A ring of symmetric functions in countably many variables is denoted by A 
(a precise definition of A see in [H). 

Let us remind some of well-known symmetric functions. 

For each m > 1 the mth power sum is 

i 

The collection of variables pi,P2, ... we denote by p. 
If we define 

Px=PxiPx2--- (2.1.5) 

for each partition A = (Ai, A2, . . .), then the px form a basic of symmetric polynomial functions 
with rational coefficients. 

Having in mind the links with soliton theory it is more convenient to consider the following 
" power sums" : 

7m=^j:^T (2.1.6) 

which is the same as Hirota-Miwa change of variables ( [1.0.81 ). The vector 

7 = (71,72,...) (2.1.7) 

is an analog of the higher times variables in KP theory. 

Complete symmetric functions are given by the generating function 

X{{l-x,z)-^ = Y.K{^)z^ (2.1.8) 

i>l n>l 

If we define h\ = h\^h\,^ ■ ■ ■ for any A, then h\ form a Z-basis of A. 

The complete symmetric functions are expressed in terms of power sums with the help of 
the following generating function: 

00 

exp ^^z^ = J2 hn{l)z'' (2.1.9) 

m=l n>l 

Now suppose that n, the number of variables finite. 
Given partition A the Schur function s\ is the quotient 



det (^x" -'J 

where denominator is the Vandermond determinant 



det fx^+"-^■^ 

^a(x) = , \ (2.1.10) 



l<i j<n 



^^K<"Oi<„,<n= n (^^-^.) (2-1-11) 



l<i,j<n 

l<i<j<n 



Both the numerator and denominator are skew-symmetric, therefore Schur function in variables 
Xi, . . . ,Xn is symmetric. 

For zero partition one puts so(x) = 1. 

The Schur functions sx{xi, . . . ,x„), where /(A) < n, form a Z-basis of An- 
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Each Schur function sx can be expressed as a polynomial in the complete symmetric func- 
tions hn- 

= det (2.1.12) 

where n > /(A). 

Orthogonality. One may define a scalar product on A by requiring that for any pair of 
partitions A and /i we have 

< Px,Pf, >= Sxf^zx (2.1.13) 

where 6xfi is a Kronecker symbol and 

where rrii = mj(A) is the number of parts of A equal to i. 
For any pair of partitions A and fi we also have 

< sx,s^>=6x,f, (2.1.15) 

so that the sx form an orthonormal basis of A. The sx such that |A| = n form an orhtonormal 
basis for all symmetric polynomials of degree n. 

2.2 A deformation of the standard scalar product and series of hy- 
pergeometric type 

Scalar product <, >r,n- Let us consider a function r which depends on a single variable n, 
the n is integer. Given partition A let us define 

rx{x) = n r{x+j-i) (2.2.1) 

Thus rx{n) is a product of functions r over all nodes of Young diagram of the partition A where 
argument of r is defined by entries i, j of a node. The value of j —i is zero on the main diagonal; 
the value j — i is called the content of the node. For instance, for the partition (3, 3, 1) the 
diagram is 



so our rx{x) is equal to r{x + 2){r{x + l))^(r(x))^r(a: — l)r{x — 2). 
For zero partition one puts tq = 1. 

Given integer n and a function on the lattice r{n),n G Z, define a scalar product, parame- 
terized by r, n as follows 

< Sx, Sf, >r,n= rx{n)5xf, (2.2.2) 

where 5x^ is Kronecker symbol. 
Let rii E Z are zeroes of r, and 

fc = min|n — nj|. (2.2.3) 
The product ( ^^2^21) is non-degenerated on A^. Really, if k = n — rii > due to the definition 



( 2.2.1 ) the factor rx{n) never vanish for the partitions which length is no more then k. In this 
case the Schur functions of k variables {sa(x'^), Z(A) < k} form a basis on A^.. li n — rii = —k < 
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then the factor rx{n) never vanish for the partitions {A : /(A') < k} , where A' is the conjugated 
partition, then {sa(x'^), /(A') < k} form a basis on A^. If r is non-vanishing function then the 
scalar product is non-degenerated on Aqo- 

Example The case r(k) = k plays a role in applications to two- matrix models. In this 
case we have a realization of this product similar to the realization of standard scalar product, 
see ( |3.3.2| ) below. We shall prove that for the symmetric functions of n variables xi, . . . , x„ we 
have the following realization 



< /(x),^(x) 



:A{d)f{d) ■ A(x)^(x) 



|x=0! 



A{d) 



n 

l<i<j<n 



d 



d 



^ dxi dxj 



(2.2.4) 



all Xi are replaced by -rr- and where r{k) = k. The proof follows from definition of the Schur 



where A(x) is the Vandermond determinant ( |2.1.11j ), f{d) is the symmetric function where 

d 

dxi 

functions ( p.l.lOD and the formula 
\ det f 9^^+"-^ ■ det f T " 



i,j=l,...,n 



|x=o = (^A/x n {K + n-i)\=<sx,s^>r,, 



l<i<j<n 

It is interesting to compare the scalar product ( |2.2.4| ) with that considered in p^.D 



(2.2.5) 



Let us generalize the example ( 2.2.4 ) 

A representation of the scalar product if r(0) = 0. Let us take 



1 d 
X = —r[D), D = X— 
X ax 



Then 



x'"■x"U=o = ^„^nr(l)r(2)■■■r(n) 
Proposition 1 We have the following representation of scalar product 



< /(x),^(x) >r,n= -rA(x)/(x) ■ A(x)^(x) |x=0, A(x) 



n 



(2.2.6) 
(2.2.7) 

(2.2.8) 



l<j<j<n 



The proof follows from 



— det (xj'^" . -det ^) , , |x=o = 5\t,rx{n) =< sx, >r,n 

77,1 \ /«j=l,...,n V / 2j=l,...,n 

A useful formula. First let us remind the so-called Caughy-Littlewood formulae 

11(1 - XiVj)-^ = SA(x)sA(y) 



(2.2.9) 



(2.2.10) 



where the number of x,y variables is irrelevant, and the Schur functions are defined via the 
formula (gXTO|) . 

If 7 = (71,72, . . .) and 7* = (71,72, • • •) be two independent semi-infinite set of variables. 
Then we have the following useful version of Cauchy-Littlewood formula: 



exp ^ m7^7;; = s xi^r) s x{l* 



(2.2.11) 



m=l 



where the sum is going over all partitions including zero partition and sx are constructed via 
(|2.1.12| ) and (|2.1.9|) . Since I did not meet the relation (|2.2.11| ) in the literature I need to 
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prove it (see ( p. 2. ID below). This proof follows from bosonic (l.li.s) and from fermionic (r.li.s) 
representations of the vacuum TL function. 

Series of hypergeometric type. Let us consider the following pair of functions of vari- 
ables Pm 

oo oo 

exp m-fmtm, exp ^ m-fmt*^ (2.2.12) 

m=l m=l 

where tm and are formal parameters. Considering pm as power sums and using formulae 
(|2.2.11| ) and ( p.2.2|) , we evaluate the scalar product of functions (|2.2.12| ): 

oo oo 

< exp m'Jrntm,exp Y ^Irnfm >r,n= Y r x{n) Sx{t) S x{t*) = ^ ^ x{n) Sx{t) S x{t*) 

m=l m=l l{X)<k A 

(2.2.13) 

where sum is going over all partitions A of the length /(A) < k, see ( 2.2.3| ), and including zero. 
The Schur functions sx(t), sxi't*) are defined with the help of 



SA(t) = det /iA,-i+i(t)i<,j</(A), exp ^ = ^ z''hk{t) (2.2.14) 

m=l k=l 

We call the series (|2.2.13| ) the series of hypergeometric type. It is not necessarily convergent 
series. If one chooses r as rational function he obtains so-called hypergeometric functions of 
matrix argument, or, in other words, hypergeometric functions related to GL{N,C)/U{N) 
zonal spherical functions [^,|^. If one takes trigonometric function r we gets q-deformed 
versions of these functions, suggested by I.G.Macdonald and studied by S.Milne in [RB . 



In case r(0) = 0, which we shall show is important in applications to models of random 
matrices, we obtain 



Tr{n, t,t*) = eS-=oEr=i*'"^rA(£) • eS™=oEr=i*m^rA(x) 



=0, X, = -r 1 (2.2.15) 



Scalar product of series of hypergeometric type. Considering the scalar product 
of functions of 7^ and using formulae (|2.2.13|) and ( 2.2.2 ), we evaluate the scalar product of 
series of hypergeometric type ( |2.2.13| ), which are related to the different choice of the function 



r (which we denote by vi and don't miss it with the notation vx of (|2.2.1| ), labelled by a 
partition!) . We find that the scalar product is a series of hypergeometric type again: 

< rr^(/c, t,7),r,.2(m,7,t*) >r,n= r,-3(0,t,t*) (2.2.16) 



where 



^sii) = f^i{k + i)r{n + i)r2{m + i) 



(2.2.17) 



Solitons, free fermions and fermionic realization of the 
scalar product 



In this section we use the language of free fermions as it was suggested in [23|, 



3.1 Fermionic operators, Fock spaces, vacuum expectations, tau 
functions []7|] 

Free fermions. We have the algebra of free fermions the Clifford algebra A over C with 
generators ipnyi'ni^ ^ ^)) satisfying the defining relations: 

[■^m, 'ipnU = [^*m, <]+ = 0; [t/jm, <]+ = ^mn- (3.1.1) 
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Any element of W = (©mezC^/^m) © (©mezC'?/'^) will be referred as a free fermion. The 
Clifford algebra has a standard representation {Fock representation) as follows. Put Wan = 
{®m<oCipm) © i®m>oCip*J, and Wcr = i®m>oCipm) © {®m<oCip;^), and consider left (resp. 
right) A-module F = A/AWan (resp F* = Wcr A\A). These are cyclic A-modules generated 
by the vectors |0) = 1 mod AWan (resp. by (0| = 1 mod WcrA), with the properties 

^m\0) = (m<0), dO) = (m>0), (3.1.2) 
(0|^^ = (m>0), (0|C = (m<0). (3.1.3) 

Vectors (0| and |0) are refereed as left and right vacuum vectors. Fermions w G Wan eliminate 
left vacuum vector, while fermions w G Wcr eliminate right vacuum vector. 

The Fock spaces F and F* are dual ones, the pairing is defined with the help of a linear 
form (0||0) on A called the vacuum expectation value. It is given by 

(0|1|0) = 1, (0|^^C|0) = 1 m<0, (0|CV'm|0) = l m>0, (3.1.4) 

{0\^m^n\0) = {OirmrjO) = 0, (0|^„dO) = m ^ u. (3.1.5) 
and by the Wick rule, which is 

{0\Wi ■ ■■W2n+l\0) = 0, {0\Wi ■ ■■W2n\0) = ^ Sgua {0\w^(^i)W^(^2)\0) ■ ■ ■ {0\Wa(2n-l)W^(2n)\0), 

(3.1.6) 

where Wk € W, and a runs over permutations such that o"(l) < cr(2), . . . ,cr(2n — 1) < (r{2n) 
and (t(1) < a{3) • • < a(2n - 1). 

Lie algebra gl{oo). Consider infinite matrices {aij)i^j^z satisfying the condition: there 
exists an such that atj = for |« — j| > A^; these matrices are called generalized Jacobian 
matrices. The generalized Jacobian matrices form Lie algebra. Lie bracket being the usual 
commutator of matrices [A, B] = AB — BA. 

Let us consider a set of linear combinations of quadratic elements ^aij : il^i'ip* :, where the 
notation :: means the normal ordering which is defined as : ipiip* := V'i^j ~ (O|'0i'0j|O)- These 
elements together with 1 span an infinite dimensional Lie algebra gl{oo): 

E ■ ^i^j •' J2 % • ^i^j ■]=J2 ■ ^i^j ■ +^0, (3.1.7) 

Cij = '^aikbkj — ^^bikttkj, (3.1.8) 

k k 

and the last term 

Co = ^ij^ji ~ (^ijbji- (3.1.9) 

i<0,j>0 i>0,j<0 

commutes with each quadratic term. 

We see that Lie algebra of quadratic elements '■ i^ii^j '■ is different of the algebra of the 
generalized Jacobian matrices; the difference is the central extension Cq. 

Bilinear identity. Now we define the operator g which is an element of the group corre- 
sponding to the Lie algebra gl{oo): 

g = expY "-nm ■ Ipntpm ■ (3.1.10) 



Using (13.1.71) it is possible to derive the following relation 

gi^n = J2'^rnAmn9, ^J^g = g^ '^nmi^*^- (3.1.11) 
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where coefficients Anm are defined by a„m- In turn ( ^.l.llD yields 



(3.1.12) 



The last relation is very important for applications in soliton theory, and it is equivalent to the 
so-called Hirota equations. 

Higher times. Let us introduce 

oo +00 +00 

Hn= ^f^rk+n, n^O, H{t) = ^ tnH^, H*{r) = ^ t^H.^. (3.1.13) 

k=~oo n=l n=l 

Hn € (7/(00), and H{t),H*{t*) belong to (7/(00) if one restricts the number of non-vanishing 
parameters tm,^^- For iJ„ we have Heisenberg algebra commutation relations: 



[Hn, Hm] = n6, 



m+n,0- 



Let us mark that 



Hn\0) = = (0|if_„, n>0 
For future purposes we introduce the following fermions 



Using ( |3.1.1| ) and ( p.l.l3| ) one obtains 



Using 



(3.1.14) 
(3.1.15) 

(3.1.16) 

(3.1.17) 
(3.1.18) 

(3.1.19) 
(3.1.20) 



^ipn®'ip*n = res^=o^(-2) ® ip*{z) 
and formulae ( |3.1.17| ),( |3.1.18| ) result in 

Therefore if g solves (|3.1.12|) then e^(*)^e^*(**) also solves ( ^.1.12|) : 

(3.1.21) 

The KP and TL tau functions First let us define vacuum vectors labelled by the integer 



res,=oV^(^) (g) ^*{z), e^(*)(7e^*(**) ® e^(*)(7e^*(**) 



M: 



{n\ = {0\^l \n) = v&^IO), 



(3.1.22) 



^'n = V'n-i---^iV'o n>0, ^^„ = <---C2d n<0, 

K = i^0^l---^n-l n>0, K = ^.itP.2---^Pn n<0. 



(3.1.23) 



Given g, which satisfy bilinear identity (|3.1.12|) , one constructs KP and TL tau-functions 
( |1.0.7| ) as follows: 

rxp(n,t) = (n|e^(*)^|n), (3.1.24) 
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TTL{n,t,t*) = (n|e^(*)(7e^*(**V)- (3.1.25) 

If one fixes the variables n, t*, then ttl niay be recognized as tkp- 

The times t = (^1,^2, • • •) and t* = (^1,^2; • • •) are called higher Toda lattice times |2^, ^ 
(the first set t is in the same time the set of higher KP times. The first times of this set ^1,^2, ^3 
are independent variables for KP equation (|1.0.5 ), which is the first nontrivial equation in the 
KP hierarchy). 

Fermions and Schur functions. In the KP theory it is suitable to use the definition of 
the Schur function s\ in terms of higher times 

SA(t) = det{hn,-i+j{t))i<ij<r, (3.1.26) 

where hm{t) is the elementary Schur polynomial defined by the Taylor's expansion: 

+00 +00 

e«(*'^) = exp(^ tkz'') = z"/i„(t). (3.1.27) 

fc=l n=0 

Lemma 1 ^ 

For —ji < ■ ■ ■ < —jk <0<is<---<ii, s — k>0 the next formula is valid: 

{s - kle^'^'^r-j, ■ ■■r-,,As ■ ■ ■ V-nlO) = (-l)^i+-+^-'=+(fc-^)(^-^+i)/25,(t), (3.1.28) 

where the partition A = {rii, . . . , n^-fc, n^-fc+i, . . . , ris-fc+jj is defined by the pair of partitions: 

(ni, . . . , Us-k) = {ii - {s - k) + l,i2 - (s - k) + 2, . . . , is^k), (3.1.29) 
{us-k+i, ■ ■ ■ ,ns-k+n) = {is-k+i, ■ ■ • - 1, . . . , jfc - 1). (3.1.30) 

The proof is achieved by direct calculation. Here (. . . | . . .) is the Frobenius notation for a 
partition, see ( |2.1.3| ). 

3.2 Partitions and free fermions 

First let us prove the formula ( |2.2.11| ). Consider the so-called vacuum TL tau function 

(0|e^(*)e^*(**)|0) (3.2.1) 



By the Heisenberg algebra commutation relation (|3.1.14| ) and by (|3.1.15| ) it is equal to 



eE™™*-*m (3.2.2) 

On the other hand one may develop e^*^**^ in Taylor series, then use the explicit form of H* 
in terms of free fermions (see (|3.1.13| ), and then use Lemma 1 of previous subsection. He gets 

E^A(t)sA(t*) (3.2.3) 

A 

where sum is going over all partitions including zero partition □. 

Lemma 2 Given partition A = {ii, . . . ,is\ji — 1, ■ ■ ■ , js — 1) let us introduce the notation 

I A) = (-l)n+-+^s^*_^^ . . . .^.jo) (3.2.4) 

(A| = (_l)^-^+-+:'-^(0|V^;^ . ■ ■ ■ -V^-,, (3.2.5) 
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Then 

(A|/x) = 6x,^ (3.2.6) 

where S is Kronecker symbol, and 

.a(H*)|0) = |A) (3.2.7) 

(0|sa(H) = (A| (3.2.8) 
here sx is defined according to (|2.2.14| ) where instead of t = (ti,t2, ■ ■ ■) vfe have 



H* = (i/_i,%...,^,...) (3.2.9) 
2 m 

H = {H„^,...,^,...) (3.2.10) 
2 m 

□ 

The proof of ( 3.2.7] ) and of (|3.2.8|) follows from using Lemma 1 and the development ( 2.2. 11|) 



ei::.i H^t- = ^ SA(H)sA(t) (3.2.11) 

A 

Lemma 3 

sa(-A)IO) =rA(0)|A) (3.2.12) 
where sx is defined according to (|2.2.14|) where instead of t = (ti, ^2, • • •) have 

A = (Ai,4^,...A,...) (3.2.13) 
2 m 

oo 

Ak= Yl rn-kMn)r{n-l)---r{n-k + l), A; = 1,2,.... (3.2.14) 

n=— oo 

□ 

The proof follows from the following. First the components of vector A mutually commute: 
[Am, Ak] = 0. Then we apply the development 

eEZ=.^rr.t^ = Y,sx{-A)sx{t) (3.2.15) 

A 

and apply the both sides to the right vacuum vector: 

eE:.i^™*"^|0) = J2sx{t)sxi-A)\0) (3.2.16) 

A 

^^^^^^ ..^ ^^.^^^^^ ^ m=i ""lO) using Taylor expansion of the exponential function and also 
use the explicit form of ([3.2.14]) , Lemma 1, the definition of |A) and also ( ]3.2.6] ). 



3.3 Vacuum expectation value as a scalar product. Symmetric func- 
tion theory consideration. 

Let us remember the definition of scalar product ( ]2.1.13] ) on the space of symmetric functions. 
It is known that one may present it directly in terms of the variables 7 and the derivatives with 
respect to these variables: 

1 1 

d = (^71' 2^12,- ■ ■ , -d^n, ■ ■ •) (3.3.1) 
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Then it is known, that the scalar product of polynomial functions, say, / and g, may be 
defined as 

{f,9) = {m-g{i)) 



1 7=0 



(3.3.2) 



For instance, using ( p.3.2|) , one easily gets 



n 



(3.3.3) 



which is a particular case of ( |2.1.13| ). 
The following statement is important 

Proposition 2 Let scalar product is defined as in ^J.^) . Then 

{f,g) = {0\f{H)g{H*)\0) 

where 



H = ( Hi, . . . , — , . . . 

n 



2 



n 



(3.3.4) 



(3.3.5) 



It follows from the fact that higher times and derivatives with respect to higher times gives 
a realization of Heisenberg algebra H^Hm — HmH^ = k6k+mfl'- 



and from the comparison of 



with 



d^^ -1 = 0, Free term of d^^ is equal to 
Hm\n) = 0, m > 0, {n\Hm = 0, m < 



(3.3.6) 

(3.3.7) 
(3.3.8) 



respectively. 

Now let us consider deformed scalar product (|2.2.2|) 



(3.3.9) 



Each polynomial function is a linear combination of Schur functions. We have the following 
realization of the deformed scalar product 



Proposition 3 

where 



{f,g)r,n={n\f{H)g{-A)\n) 



t± — -Hi, — , • • • , , 

2 m 



A= (Ai,^,...,^. 



(3.3.10) 
(3.3.11) 



Hk and are defined respectively by ( \3.1.13{ ) and (\3.2.1^ . 
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3.4 KP tau-function rr{n,t,V) 

For the collection of independent variables t* = (^1,^2; ■ ■ ■) '^Q denote 

00 

m=l 



(3.4.1) 



where were defined above in ( p.2.14| ). 

Using the explicit form of Am ( p.2.14|) and anti-commutation relations ( p.l.l|) we obtain 



where operators 



+00 .1 

m=l 



D = z^, r'{D) = r{-D) 



(3.4.2) 
(3.4.3) 

(3.4.4) 



act on all functions of z on the right hand side according to the rule which was described in the 
beginning of the subsection 3.4, namely r{D)-z^ = r{n)z^. The exponents in ( |3.4.2| ),( f3.4.3| ),( P".4.4| ) 
are considered as their Taylor series. 

Using relations (|3.4.2|) , (|3.4.3| ) and the fact that inside res^ the operator -r{D) is the con- 
jugate of ir(— D) = jr'{D), we get 

Lemma The fermionic operator e^'**^ solves the bilinear identity ( p.l.l2| ): 



ie'Sz=oip{z) ® ip*{z), e 



(3.4.5) 



□ 



Of cause this Lemma is also a result of a general treating in [0. 
Now it is natural to consider the following tau function: 

r,(n, t,t*) := (nje^^^^e-^^** V)- (3-4.6) 
Let us use the fermionic realization of scalar product ( p. 3. 10 ). We immediately obtain 



Proposition 4 



r,(n,t,r) = (eS^=i™*'"^'",eS™=i'"*-^'"),,„ 



(3.4.7) 



Due to ( p.2.131) we get 



Proposition 5 



rr(n, t,t*) = Yrx{n)sx{t)sx{t'' 

A 



(3.4.8) 



We shall not consider the problem of convergence of this series. The variables M, t play 
the role of KP higher times, t* is a collection of group times for a commuting subalgebra of 
additional symmetries of KP (see |3^, |13| and Remark 7 in [0). From different point of 
view ( |3.4.(^ ) is a tau-function of two-dimensional Toda lattice with two sets of continuous 
variables t, t* and one discrete variable M. Formula (|3.4.8| ) is symmetric with respect to t ^ t*. 
For given r we define the function r': 



r'[n) := r{—n). 



(3.4.9) 
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Tau function r,. has properties 



Tr{n, t, t*) = Tr{n, t*, t), Tr'{—n, — t, — t*) = Tr{n, t, t*). 



(3.4.10) 



Also Tr{n, t, t*) does not change if tm oJ^tm, ^ cl "^t^, m = 1, 2, . . . ,, and it does change 



if t 



a'^tm,m = 1,2, . . .,r{n) 



a ^r{n). ( |3.4.1(]| ) follows from the relations 
r'xin) = ry{-n), Sx{t) = {-f^sy{-t). (3.4.11) 
Let us notice that tau-function ( ^.4.8] ) can be viewed as a result of action of additional 



symmetries |3^,p0|,||3^ on the vacuum tau-function . 



3.5 Algebra of "^DO on a circle and a development of e"^^* \ e"^^*^ 
Given functions r, f, the operators 



- J2 rin)■■■r{n-m + l)^|Jn^p^_^, m=l,2, 



(3.5.1) 



m 



1,2, 



(3.5.2) 



belong to the Lie algebra of pseudo-differential operators on a circle with central extension. 
These operators may be also rewritten in the following form 



Ar, 



271 



i=^:^*(^) Qr(Z}))"-^(z) 



(3.5.3) 



(3.5.4) 



where D 



pseudo-differential operators r{D),r{D) act on functions of z to the right. 



This action is given as follows: r(D) ■ z^ = r{n)z'^^ where z^,n E Z is the basis of holomorphic 
functions in the punctured disk 1 > \z\ > 0. Central extensions of the Lie algebra of generators 
( 3.5.3| ) (remember the subsection 3.1) are described by the formulae 



Un{Am, Ak) = Smkr{n + m - 1) ■ ■ ■ r{n)r{n) ■ ■ ■ r{n - m + 1), Un - uj„ 







(3.5.5) 



where n is an integer. These extension which differs by the choice of n are cohomological 
ones. The choice corresponds to the choice of normal ordering ::, which may chosen in a 
different way as : a : := a — {n\a\n), n is an integer. (We choose it by n = see subsection 
3.1). When functions r, f are polynomials the operators Qr(D)^ , (f(D)z)™' belong to the 
Woo algebra, while the fermionic operators (|3.2.14|) ,( p.5.2|) belong to the algebra with central 
extension denoted by Woo Pi • 



We are interested in calculating the vacuum expectation value of different products of op- 
erators of type e^"^"'*'" and e^"^'"*'". 

Lemma ||12|| Let partitions A = (zi, . . . , is\ji — l, ■ ■ ■ , js — 1) and = (ii, . . . , — 1, • • • ,jr— 
1) satisfy the relation A 3 yU. The following is valid: 



mi ■ ■ ■ ■ ■ ■ ^_j,e^(** V-,, ■ ■ ■ ■■■^n |0) = 

= (-l)^'^+-+^''-+^'^+-+^'^SA/^(t*)rA/;.(0) 



(3.5.6) 
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(o|^;^ ■ ■ ■ ■ ■ ■ ^-,,e^(*V!j, ■ ■ ■ ^!,v^i. ■ ■ ■ ^i, |o) = 

= (_l)ii+-+>+.-i+-+>sv^(t)f,/^(0) 
where sx/^{t) is a skew Schur function 

SA//.(t) = det (/iA.„^^_i+j(t) 

and 



;n) := [| r(ra + j-ii), rA/^(n) := r{n + j-i) 



(3.5.7) 

(3.5.8) 
(3.5.9) 



□ 



The proof is achieved by a development ofe'^ = l + A + -- -,e^ = l + A + -- - and the di- 
rect evaluation of vacuum expectations ( p.5.6| ),( p.5.7| ), and a determinant formula for the skew 
Schur function of Example 22 in Sec 5 of f^. 

Let us introduce vectors 

\X,n) = (-l)^l+---+^'=C,i+n---C,.+n^i.+n---V^n+n|ri) (3.5.10) 

{\,n\ = (-l)^'^+---+^-^(n|^;^+„---V^;^+„^_,,+„---^_,,+„ (3.5.11) 

As we see 

{\,n\fx,m) = SmnhfM (3.5.12) 

We see that Fock vectors (|3.5.1CI|) (and (|3.5.11 )) form a ortho normal basis in the Fock spaces 
F (and respectively) F*. 

Remark. From (|3.1.16| ) we obtain 



f-- iijizi)---ij{zn)\0><0\r{zn)---r{zi)= E |A,n><A,n| (3.5.13) 

X,l{X)<n 

which projects any state to the component F„ of the Fock space F = ^nez Fn ■ ^ 
Using Lemma 1 and the development ( |3.2.15| ) we get 
Lemma 

|A,n) = SA(H*)|n), (A, n| = (n|sA(H) 

and 

Sx{-A)\n) = rA(n)sA(H*)|n), {n\sx{k) = fx{,n){n\sxili) 



Hr, 



2 m 



~ ( 7 A2 Am 

\ 2 m 



(3.5.14) 

(3.5.15) 
(3.5.16) 



where Am, Am are given by ( |3.2.14| ), (|3X^ ). □ 
There are the following developments: 



Proposition 6 



-A(t*) 



E E \^^,n)sx/^,{t*)rx/^,{n){X, 



n\ 



neZ XDfM 



n\ 



(3.5.17) 
(3.5.18) 



nSZ XD^ 

where SA/^(t) is the skew Schur function ^3. 5. d^ ) and rx/^{n),fx/ii{n) are defined in ( \3. 5.9[ ). 
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The proof follows from the relation ^ 

(■Sa/m, s^) = {sx, s^s^) (3.5.19) 
Remark. It follows from (|3.5.17] ), ( p.5.18|) 



s,i-A) = J2Y.C;urx/,{n)\l^^^){>^^^l ^M) = E E C;.rv,.H|A, n| (3.5.20) 
where coefficients C^^ define the decomposition of the product of Schur functions s^s^ = 

Sa C'^u^X- 

3.6 Schur function expansion of tau function. Loop scalar products 

According to formulae (|3.1.25|) and ( p.5.12| ) we have 



r,(n,t,r) = (n|e^(*)^?e^(**)|n) =Y^Sx{t)gx,{n)s,{V) (3.6.1) 

X,ij, 

where is the matrix element 

gxi^iji) = {\n\g\ii,n) (3.6.2) 

Therefore we have the property: 

{gi92)xt,{n) = Y^{gi)xv{n){g2)y^,{n) (3.6.3) 

We can introduce the following 'integrable' scalar product 

< Sx, >g,n= gxf,{n) (3.6.4) 

In general this form is a degenerated one. 
We see that we obviously have 

< T-<;i(ni,t,7),r<;j(n2,7,t*) >g,„= rg3(0,t,t*), gs{k) = gi{ni + k)g{n + k)g2{n2 + k) (3.6.5) 

In the previous consideration we actually took 

A-AmH = ^Xi,rx{n) (3.6.6) 

Different examples of non diagonal gx^i expressed in terms of skew Schur functions one can 
obtain using the Proposition of previous subsection. It will be used in the section "Further 
generalization" below. Let us note that in fact non diagonal gx^ were constructed in the form 
of certain determinants in the papers [|28|,[|^,|^ • 



Now let us consider the following scalar product of Tgi with itself: 

« VW »a,n--=< v(n,7,7) >3,„= E^Am(^)vH (3-6-7) 

X,n 

One can easily see that this scalar product is the trace of operator g in the component of the 
Fock space of the charge ra, namely, F„: 

« VW »g,n=< Tg>{n,'y,'y) >g,n= J2 9x,M9'^xin) = E(^' ^l^'^l^' ^) = TraceF„(^'^) 

X,fi X 

(3.6.8) 

Let us consider the following loop scalar product: 

(3.6.9) 

E ^lAiAaH • ■ -^PApAiH = TraceF„(^i^^2^- • - ^p^) (3.6.10) 

Ai,...,Ap 

As we shall see below the loop scalar product describes certain matrix integrals. 
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4 Angle integration of tau functions of matrix argument. 
Integration over complex matrices 

In this section we shall generalize some known solvable matrix integrals. By the solvable matrix 
integrals we mean those which can be presented as tau functions and which can be evaluated by 
the method of orthogonal polynomials. Let us make a reference to the related paper where 
by solvable matrix model the author mean the model which admits a reduction of number of 
integrals (from the order ~A^^ to the order A^) via a character expansion method. Let us refer 



to the papers [35|,|48],[34 where different generalizations of HCIZ integral were considered. 



4.1 Useful formulae 

In this section we shall exploit the following formulae of an integration of Schur functions over 
the unitary group 0]. First is the formula 

/ sMUBU-^)d.U^'-^^^^, (4.1.1) 

Ju{n) Sx{In) 

where d^:U is the Haar measure on the U{n). Then we have 

/ sx{AU)s,{U-'B)dM = ^^'^Am- (4.1.2) 

Ju(n) Sx{In) 

If we have complex matrices Z then there are formulae B 



TT 



sMZBZ^)e-^^^'' n d^Z,,d-sZ,, = "'^f'i^^ , (4.1.3) 
J-oo - j^^ SaI^I, u, U, U, . . .) 



and 



TT 



..(AZ).,(Z+5)e^Trzz+ -q ^sr^^^.^cj^^^. ^ /^>^(AB) ^^^^^^^ 

J-oo - -^^ SA(,i, U, U, U, . . .J 



Then we need 



sa(1, 0, 0,0,...) = Hx = + a; - z - J + 1) (4.1.5) 



where H\ is the hook product. The Pochhammer's symbol related to a partition A = (Ai, . . . , A^) 
is the following product of the Pochhammer's symbols (a) a := {a)x^{a — 1)a2 ■ ■ ■ {a — k + ^)\^- 
Then 

(a)A = i^ASA(t(a)), Ha = H^s^^h) (4.1.6) 

In addition we have the relation (bosonic (l.h.s) and fermionic (r.h.s) representation of the 
vacuum TL tau function) 

oo 

exp mt^t*m = E ^A(t)sA(t*), (4.1.7) 

m=l X 

which is a general version of Cauchy-Littlewood identity 
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Let t:, = (l,0,0,0,...),t(a) = (f , |, |, . . .), r(fc) = nLi(A; + a.) ni=i(A; + /(A) is the 

number of nonvanishing parts of A. Let us write down the formulae for the hypergeometric 



function of a matrix argument i2t]^ see Appendix C: 

/ai + M, ...,ap + M 



p-^ s 



bi + M,---,bs + M 



X\ = r,(M, t(x"),t: 



A 

i(A)<n 



nu sx{t{h + M)) 



(4.1.8) 



Now let r{k) = HiLil^ + ai){k + n — M) ^Y\i^i{k + bi) ^. The hypergeometric functions of 
two matrix arguments j^/, see Appendix C, is 



'ai + M, . . . ,ap + M 
61 + M, ■ ■ ■ , 6, + M 



E 771?]— ^r7Tv(*A(l, 0,0,0,...)) 



X,rj =r,(M,t(x"),r(y")) 
p+l SA(x")sA(y' 



A 

!(A)<n 



nLi^A(t(fefe + M)) 



(4.1.9) 



Here x" = (xi, . . . , a;„), y" = (?/i, . . . , ?/„) are eigenvalues of matrices X, "K. 

Let us mark also the determinant representations of (|4.1.8|) and of ( [4.1.9|) (details 
we find in Appendix A), which are of importance in applications of the method of orthogonal 
polynomials to the new matrix models, which we shall consider below. 

The determinant representations have forms (the reader may find details in p| , [pJ]] ) : 



r^(M, t(x"),t* 



and 



where 



r,(M,t(x"),t*(y")) 



det (xrV,(M- A; + l,t(xi),t*))"^_^ 
A(x") 

cj^ det(r^(M-n + l,Xi,|/j))". ^ 



A(x"l 



A(x")A(y") 
det(xr^). 



i<3 



n-1 



Cn = n (K^)) 



k—n 



(4.1.10) 

(4.1.11) 
(4.1.12) 
(4.1.13) 



k=0 



Tr{M -n + 1, Xi, Vj) = 1 + r(M -n + l)xiyj + r{M -n + l)r(M -n + 2)x-y- + ■■■ (4.1.14) 
Below we need mainly the case M = n. 



4.2 Certain matrix integrals and expansion in Schur functions 

In this subsection we shall consider few integrals of exponential functions. We shall pick up the 
solvable models. 

Integration over complex matrices. Different integrals over complex matrices were 



considered in the papers of 1. Kostov, see example. 

(A) Using relations (|4.1.7|) and ( |4.1.3|) we evaluate the following integral over a complex 
matrix M as the expansion over Schur functions 



sx{t)sx{A)sx{B) 



x,iw<n ^a(1, 0,0,0, 



(4.2.1) 
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where (fZ = Tr"" Ulj=i dMZijd'^Zij. 

Let us compare this series with the relation ( |4.1.9| ). To identify them we have the follow- 



ing prescription. He will specify at least one of the arguments of the Schur functions in the 
numerator with one of the sets: with (1, 0, 0, 0, . . .) or with (j, |, |, . . .). Let us notice that 
if we take a = n we get the Schur function of the unit matrix /„, which is the dimension of 
the representation labeled by A. Also we mark that sa(1, 0, 0, 0, . . .) is the inverse of the hook 
product Hx- 

Thus we have few opportunities. 

(Al) First, putting t = (1, 0, 0,0,.. .) we get 



|eS:^.Tr(AZBZ.)^_TrzzVz= sM)s.{B) =fl{l - a^h,)'' = ,:F,{n\A, B) 

X,l{X)<n i,j 

(4.2.2) 

(A2) Second, putting t(a) = {j, f , f , • • < 1, we get the expansion of the following 

matrix integral 

[ det - ^AZBzA e-^^^^^£Z = J2 ia)xSxiA)sx{B) = (a, n\A, B) (4.2.3) 

which in general is a divergent series. 

(A3) Then let us choose the matrix B as diagonal matrix with 5fcfc = = l,...,n. Let 
us consider the limit n ^ oo. In the limit with the help of ( |4.1.5|) we get the hypergeometric 
function of matrix argument 

E («)A^ = iFo(a;A) (4.2.4) 

\,l{\)<n 



(B) With the help of ( |4.L7| ), ( |4.L4| ) we get the development 



= / eYZ^. t„Tr(AZ)-+E:.i cTr(z+B)'" -Trzz+^2 z= ^ sx{t)sx(X)sx{AB) ^ 
J sa(1, 0,0,0,...) ^ • • ^ 



\,l{X)<n 



Comparing it with (^4.1.9|) we obtain three solvable cases. 



(Bl) The first one is t* = (1, 0, 0, 0, . . .). Then we have 

/^ = ei::^i*-Tr(ABr (4.2.6) 

(B2) The second one is AB = In, then we get 

Ib= (n)ASA(t)sA(t*) =2^o(n,n|t,r) (4.2.7) 

\,l{X)<n 

This expansion coincides with the expansion for the model of two Hermitian random matrices 
[?]• 

(B3) The third one is t* = (f , f , |, . . .), then 

/b= E ia)xSxit)sx{AB) =2J'o{a,n\t,AB) (4.2.8) 

X,l{X)<n 

Integration over unitary matrices. 
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(C) The combination of the formula ( |4.1.1D with the formula ( [4.1.7D yields Schur function 
representation for different matrix integral: 

/ e^^-^^-^'i^^^^TdM= E '^^'^'^1^1'^^^^ (4.2.9) 



Similar series were considered in the papers [55|,[5G]. 

Looking at this formula we see that one can choose t in such a way, that he gets the 
hypergeometric functions of matrix arguments (|4.1.9| ). 

As a consequence we evaluate the following matrix integrals 

(CI) HCIZ integral 

TT(AUBU+),rr X- ■^xiA)sxiB) ^ I . r.^ det(e"'^0 o 

e )d^U= 2^ — — — = oJ^o [A,B)=Cn (4.2.10) 

u(n) Hxsxiln) A(a)A(6) 

(which we obtain putting ti = = 0,m > 1) and which is the hypergeometric function of 
matrix arguments (|4.1.9| ). For the last equality see ( [4.1.11|) . 
(C2) 



•JuH V ^ ^{^^^ HxSxiln) A(a)A(6) 

(4.2.11) 

which we obtain by putting tm = —a/m, m = 1, 2, . . .. For the last equality see (|4.1.11|) . 

(D) Let us remember the following integral , which is used for the study of two-dimensional 
QCD ||68| and called Gross- Witten one plaquette model 

= / e-^T^(^+^^)rf,f/ (4.2.12) 

Juin) 



\n—l—a 



About this integral see also |59 



A generalization of this integral was considered in |^ : 

/G^(J,J+)= / eTr(^^+^^^+)4f/ (4.2.13) 

JU{n) 

where the method of orthogonal polynomials was applied and the links with the KP equation 
and with the so-called generalized Kontsevich model was shown. 

Using the orthogonality of Schur functions we easily get the development 

/^^(J,J+)= E '^^^'^'^^y.^'°'°'---^ =o^o(jJ+;(l,0,0,0,...)) (4.2.14) 

\,l(\)<n v"'''^ 

which is the TL tau function of hypergeometric type. 

Let us further generalize the matrix integral ( ^.2.13| ) as follows 



/ ei:..*-Tr(A^r+E.cTr(^-Br^^^^ .,(Ag).,(t).,(r) ^ 15) 

Now we consider two special cases of the last integral. 

(Dl) For the first example we take t*^ = h/m,m = 1, 2, . . . < 1), then the integral ([4.2. 15|) 

is 

/ ei:™*-Tr(A^)'"det(l-f/-i5)-'4f/= ^ (^) /A(Ai^)3A(t) ^ 

(4.2.16) 
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By ( [4.1. 101 ) we have the determinant representation of ( [4.2.161 ). 

(D2) For the second example we take tm = a/m,tj^ = b/m,m 
the integral ( 4.2.15 ) is 

det (1 - AUy^ det (l - U-^b) d,U : 



1, 2, . . . (3?a, 3^6 < 1), then 



U{n) 



E 

X,l{\)<n 



{a)x{h)xsMB) 
(n)A 



= 2^1 {aMn\AB) 

(4.2.17) 

which is Gauss hypergeometric function of the matrix argument AB. By ( [4.1.10[) we have 
the determinant representation of ( 4.2.17[ ). Let us mark that for integer a, 6 the logarithmic 
derivative of the hypergeometric function 2-^1 {a,h\n\AB) solves Painleve V equation |T^. 

(D3) The third example is well-known model of unitary matrices {A 
about this model. Now we present the following representation 



B = 1). See 



U{n) 



X,l{X)<n 



(4.2.18) 



Due to the restriction Z(A) < n it is not the l.h.s of Cauchy-Littlewood formula ([4.1.7[ ). Let 
us mark that the right hand side of ( [4.2.18[ ) is the subject of the so-called Gessel theorem if 
one take t = t(x™), t* = t*(y™), n < m. This case was considered in [15|. In this case we 
have different determinant representations by ( |4.1.10D , by ( [4.1.11[ ) and also by the formula for 
semi-infinite TL tau-function (|6.8.1|) . 



4.3 Tau functions of a matrix argument and angle integration. The 
integration of tau functions over complex matrices 

It is suitable to introduce the notion of tau function of matrix argument by analogy with the 
hypergeometric function of matrix argument ( 4.1.8[) ,( ^.1.9[) : 



r (n, X, r ) := r (n, t(x"), t* ) , r (n, X, Y) := r (n, t(x"), t*(y")) (4.3.1) 

where Hirota-Miwa variables xi, . . . , are just eigenvalues of a matrix X. We shall use large 
letters for this matrix argument. As it is in the case of the hypergeometric functions the tau 
function by definition depends only on eigenvalues of the matrix. 
Generalization of HCIZ integral 

Let us consider the averaging of the TL tau function Tr{n,XUYU^ ,t*): 

I={Tr(n,XUYU\t*))u:= [ Tr (n, XUYU\t*) dM = E ''^^^'j^\''^^*K ,in) 

(4.3.2) 

If we take t* = (^j, ^, . . we get the following TL tau function Tf{n, X, Y) 

(4.3.3) 

which is certain hypergeometric function ( [1.1. 9[ ) in case r rational. The r.h.s is the determinant 
representation of the l.h.s. , see ( [1.1.11[) . A(x) is the Vandermond determinant. 
Generalization of the integral over complex matrices 

Let us consider the averaging of the TL tau function Tr{n, AZBZ^ ,t*): 

Ic = (r, (n,XZYZ+,t*))z := J r, (n,XZYZ+,t*) e-^^^^*d^Z (4.3.4) 
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We get 



lr.{n,XZYZ\t*)e-^^''U'Z= ^ ^^^^^^^^-aH (4-3.5) 



X,l(X)<n 



sa 1,0,0, 



If we take t* = i^j, ^, . . .j we get new TL tau function Tf{n, X, Y), where 

f(m) = (a + m — l)r(m) 

Finely we obtain 

Ic = ^f('^,x,y) = 



det{Tf{l,Xi,yj)}\ 



i,j=l,...,n 



A{x)A{y) 
see (^TTT|) . 

Generalization of Gross- Witten integral. This is the integral 



(4.3.6) 
(4.3.7) 



I{XY)= [ Tr{n,t,XU)Tf(n,U~^YX)d,U= ^ 



sx{XY)sx{t)s^{V 



\,l{X)<n 



-rx[n)rx[n) 

(4.3.8) 

For rr = 1 the integral ( |4.3.8| ) is equal to the integral ([4.2. 15|) . 

Let Zi,i = 1, . . . ,n are eigenvalues of the matrix XY. The matrix integral ( |4.3.8D is TL tau 
function if 

(a) = a/ m,m = 1,2, .. . (a is any complex number), then 



I{XY)= J2 {a)xrx{n)rx{n 

\,l{\)<n 



sx{XY)sx{t) _ det {zr\'{n - k + l,t{z,),t*)) H,^, 



Ha 



Afz") 



where r'(m) = '^'^^_\ r{m)r{m). The last equality is due to ( |4.1.10| ). 

(b) t*^ = 6i^rn, then 

Tfv^^ ^ sx{XY)sx{t) . ^ . det(zrV(rz-fc + l,t(z,),t*; 

I{XY) = 2^ — rx{n)rx{n) - 



(4.3.9) 



n 

U,fc=l 



X,l{X)<n 



Ha 



-r[m}r[m 



A(z") 

. The last equality is due to ( |4.1.10| ). 



(4.3.10) 



where r'(m) = — , 

(c) and (d) one gets by t t*. 
(e) XY = In then 

I{XY) = Trf{n,t,t*) 

This answer may be also derived from ( |3.5.13| ). 

Remark. Let us notice that the following averaging over unitary group U{n) 

sx{XY) 



(4.3.11) 



I{XY) = I Tr (n,XU,U-^Y) dM = VtaI 

Ju{n) ^ ^ V 



n] 



Sx{In) 



(4.3.12) 



is not TL tau function. □ 

In the similar way one can consider the integral over complex matrix Z 



sxiXY)sxit)sx{t*) ~ 



I{XY) = frr [n, t, XZ) r, {n, Z+Y, t*) e-^'^^^U'Z = ^ I^^:^^Zr^Aip^f,(n)r,(n) 

x,i{x)<n -^aU, u, U, U, . . .J 

(4.3.13) 



The reader will pick up the solvable cases in the way similar to the previous cases. 
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5 Integral representation of the scalar product and the 
integrals over eigenvalues of random matrices 

5.1 Integral representation of the scalar product <, >r^n 

In case the function r has zero the scalar product (|9.1.1|) is degenerate one. For simphcity let 
us take r(0) = 0, r{k) 0, k < n. The scalar product ( |9.1.1j ) non generate on the subspace of 
symmetric functions spanned by Schur functions {s\,l{X) < n}, /(A) is the length of partition 
A. 

To find the integral representation of this scalar product it is enough to find a function /i of 
one variable with the following property. If it depends on the product of two variables x and 
y, then moments of this function have a diagonal form: 

where F is some integration domain to be defined later, and values of r{k) are certain given 
numbers, and A is finite and does not vanish: 

A^oo, A^O, (5.1.2) 

The appropriate /i we shall denote by /i^. 

Example 1. The variable y is a complex conjugated of x = z. The integration domain is 
the whole complex plane, we suggest that ( ^.1.2 ) is fulfilled. Then obviously 



z''z"'fXr{\z\'^)dzdz = ASnm.ril)r{2)---r{n), n,m = 1,2,... (5.1.3) 

Example 2. The variable x is real, and the variable y is imaginary one. The integration 
is going over the real and over the imaginary axes. We take Urixy) = e~^^ and obtain that 
r{n) = n: 

' x'^y'^e-^'ydxdy = -2n^/^5n„^n\, n, m = 0, 1, 2, . . . (5.1.4) 



Example 3. Both variables x, y lie on the unit circle: x = e^'^'^, y = e^'^'^ . We take 
\ir{xy) = {xy)~^e^^'^^ and obtain that r(n) = 1/n: 

x-y^e^-y^-"^ = -AT^^r^m-,. n,m = 0,1,2,... (5.1.5) 
xy n\ 

□ 



Remark. We can find an appropriate function /i,. in a following way. Let r be an analytical 
function of one variable and consider r{Dx), = x-^ which acts on the power functions x'^ by 
the multiplication by r(s). One suggests that there exists a function fir{x) with the following 
properties 

(l-^r(-D,))/i,(x) = 0, D^ = x^ (5.1.6) 
One may find solutions by using the Mellin transformation as 

/i,(x) = / x'e-^^~'"^Ms (5.1.7) 
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where 



(5.1.8) 

The integration domain F in (|5.1.1|) should be chosen in such a way that the operator 
^r{D), D = x-^ is conjugated to the operator ^r(—D), D = x^: 



-r{D) 



X 



-r(-D] 



X 



on the appropriate class of functions, as used in ( |5.1.1| ) below. 
At last with the help of (|5T^ ) and (|T|) we get (|5TT|) . □ 

In case r(0) = the series 



(5.1.9) 



IJ.r{x) = 1 + 



X 



+ 



X 



r -1 r -1 r -2 



+ 



(5.1.10) 



is the formal solution of ( |5.1.6| ). This series is equal to the tau-function of hypergeometric type 
related to 1/r' , where 

r'{n) := r{~n) (5.1.11) 



fir{xy) = n/r'{l,t{x),t*{y)), kh = x^, ktl = y^ (5.1.12) 

Let us suggest that the series ( |5.1.10| ) is convergent (this fact depends on the choice of the 
function r) , and the integral ( |5.1.1|) exists {A ^ 0, oo), then one may choose this /ij. for 
(|X23). 

Example la. We consider the case of Example 1 above. The hypergeometric function 
pFs(ai, . . . , 6i, . . . , hg] —\z\'^)^p < s may be taken as an example of the function /i,, of the 
Example 1. We put |zp = x below. We have 

r x%F,{a,, ...,a„h,...,bs; -x)dx = n^=i(l ' ^On+i ^ ^g_^_^3^ 
see for instance |0. In this case 

It is known that the hypergeometric function pFg ,p < s 

{a„ ...,a„h,..., b,; -x) = £ ^^^^^^ ' ' ^^P^-f , (5.1.15) 

is a convergent hypergeometric series. It is of form ( p.l.lO| ). This series is a solution of 

(x + D^^^iy^^V"^ TT 1 p^s{ai, . . . , ap; 61, . . . , 6,; -x) = (5.1.16) 

which the hypergeometric equation ( |8.0.12| ) rewritten in the form ( [5.1.6| ), where xy is replaced 
by —X. The operator r{—Dx) acts on each a;" of the convergent series as the multiplication by 
r(— n).n 
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Example lb. We choose r{k) = —k{a— 1 — k)^^, then iFo{a; x) = (1 — x)""" is the solution 
of ( p.l.6|) , for the case 3?a < the condition (|5.1.9| ) is fulfilled inside the interval [0, 1], and for 
< 1 we have 



x"i-Fo(a; x)dx = A 



{2 -a) 



A 



(5.1.17) 



□ 



Now we shall turn to the multiply integrals. Let us evaluate the integral 



A- 



A(x)A(y)sA(x)sj,(y) W fJ^r{xkyk)dxkdyk 



k=l 



where x = (xi, . . . ,x„),y = {yi, ...,yn) 



A(z) 



HiZi - Zj) 



(5.1.18) 



(5.1.19) 



Using ( 5.1.1|) , the definition of the Schur functions ( 2.1.10 ) and the definition of rx{n) ( 2.2.1|) 
(see PI for details) we finely get the relation 



(5.1.20) 



and we obtain 
Proposition 7 

{f,9)r,n = ^ 

Then it follows that 
Proposition 8 

{f,9)r,n 



A(x)A(y)/(x)5((y) f^rix kyk)dxkdyk 



(5.1.21) 



k=l 



1 



/ ■■■ / /(x)(7(y)det(Mx.l/,)),,=i,..„„A(x)A(y)nrf^A (5.1.22) 

The proof of the last Proposition follows from the fact that the function det {fi{xiyj))- ^ 
is an anti-symmetric function in variables Xi, . . . Then by changing notations of the vari- 
ables Xi inside of the integral we obtain the Proposition. 

Let us remember that if 

n n 

t(x) = (ti(x),t2(x),...), t*(y) = (tKy),t;(y),...), mtm = j:^T, nif^ = EyT 



i=l 



4 = 1 



(5.L23) 

we have the following determinant representation , |jTO| (see an Appendix below ( |6.7.4| ), 

(iron) , 



n-l 



A(x)A(y)rr {n, t(x), t*(y)) = c„det (r^ (1, x„ yk))lk=i, = H i^i^)) 



k—n 



(5.1.24) 



k=0 



where x = (xi, . . . , x„), y = (j/i, . . . , yn) which is true on the level of formal series. 

Ti/r' (n, t(x), t*(y)) is the tau-function ( |3.4.8|) , where the function r is chosen as 1/r' 



n/r' (n,t(x),t*(y)) = 1 + (^"=i^';)(p=i^'- + 



r -1 



(5.1.25) 



Due to determinant representation ( 5.1.24 ) we see that this series is a convergent one, since we 
consider the series (|5.1.12| ) to be convergent. 
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Proposition 9 Let the tau function ri/^' of \3.4-^ , where r'{n) := r{—n), is a convergent 
series, and let r(0) = 0. Then we have the following realization of scalar product (}^.2.2[ ) 

{f,g)r,n = A--c-' j ■■■ j /(x)^7(y)rv.'(n,t(x),t*(y))(A(x))2(A(y)f ndx.dy, (5.1.26) 

-'^ -'^ k=i 

The proof follows from the fact that the function Tij^i [n, t(x), t*(y)) is a symmetric function 
in variables xi, . . . , x„. Then by changing notations of the variables Xi inside of the integral we 
obtain the Proposition. 

Then it follows from ( p.l.26| ),( |2.2.13| ) that for the same conditions the following duality 



relation is true 

Tr{n,t,t*))= (5.1.27) 

A-'^c-' f ■■■/ eEL.E:..(*™-r+CC)ri/,,(n,t(x),t*(y))(A(x))2(A(y))^nrfa:fcrfy, (5.1.28) 
•'^ •'^ k=i 

Also due to ( p.2.16|) , (p.2.17|) we have the following 

r^3(0,t,t*)= (5.1.29) 

A-X' [ ■■■ [ '^n (A;,t,t*(x))ri/,, (n,t(x),t*(y))r,, {m, t {y ), t* ) A\^)A\y) f[ dx^dy^ 
•^^ -^^ k=i 



(5.1.30) 

where 

fsi'i) = + i)r2{m + i)ri^n + i) (5.1.31) 
5.2 Normal matrix model 

Let us consider a model of normal matrices. A matrix is called normal if it commutes with 
Hermitian conjugated matrix. 

Let M be n by n normal matrix. The following integral is called the model of normal 
matrices: 

/^*^(n, t, t*; u) = J ciMdM+e^(^*^^)+^^ = (5.2.1) 

C[ ■■■[ |A(z)|'eELif^(-»-"0+^i(-0+^2(-">)fT(;^^(;^^ (5.2.2) 
Jc Jc fj^ 

where dM = ni<fc d^Mikd^Mik Dili dMa^ in the last integral the integration is going over com- 
plex planes of eigenvalues C is a number which appears due to the matrix 'angle integration' 
over V (M = VZV^^, where Z = diag{zn}), and U is defined by u = {ui,U2, ■ ■ .): 

oo oo oo 

U{MM+) = u„Tr(MM+)", Vi{M) = J2 ^™TrM"^, V2{M) = ^Tr (m+)" 

ra=l m=l m=l 

(5.2.3) 

n 

Aiz) = Hizi - Zk) (5.2.4) 

i<k 

and A{z) = 1 for n = 1. 

As we see the integral ( p.2.2| ) has a form of ( p.l.21| ), where x = z, y = z , see Example 1 of 
this section, and 

/(z) = e^tiEZ^i'rnzr ^ ^(2) = eELiE™=i*^-T, i^r{ziZi) = e^^'^'^^ (5.2.5) 



28 



Now let us compare e'^ with fir of ( p. 1.101 ). We choose the integration domain F to be the whole 
complex plane. We suppose that 



^ / eSn=i """"rf^ = A ^ 0, oo 



(5.2.6) 



The condition (|5.1.2|) is fulfilled. If he function e*^*^'^' ^ has a form of ( |5.1.10|) (each Taylor 
coefficient is non- vanishing) , then the set of variables u = {ui,U2,...) is related to the set 
{r{n),n = —1, —2, . . .} as follows 



/im fu) 



r(— 1) r{—m) 
Example. = Then 



r — m 



^ — 1- gZ^fe=i 

/im(u) 



m=0 



r(m+ 11 



r{—m) = m 



(5.2.7) 



(5.2.8) 



□ 



Finely formula ( p.2.13|) , where Pm = J27=i yields the series for (|5.2.1D : 



J^*'(n,t,t*;u) =< eS"=iS^=i*'"^«",eSr=iE:=iC^r ^ rA(n)sA(t)sA(t*) (5.2.9) 



5.3 Two-matrix model 

Details of this and of the next section one can find in 0. 

Let us evaluate the following integral over n hj n matrices Mi and M2, where Mi is a 
Hermitian matrix and M2 is an anti-Hermitian one 

P^'^'in, t, t*) = J e^^(^^^)+^^(*^^)e-Tr^^^^^ciMidM2 (5.3.1) 



Vl(Mi) = V t^TrMf, V2{M2) = V CTrM^"^ (5.3.2) 



where 

00 00 
^t^TrMf, V2{M2)=Y.^n.- 

m=l m=l 

It is well-known ||16[)@] that this integral reduces to the integral over eigenvalues Xi and yi of 
matrices Mi and M2 respectively. 

To do it first we diagonalize the matrices: Mi = UiXUi^^M2 = y/ —1U2YU2^ , where 
Ui,U2 are unitary matrices and X,Y are diagonal matrices with real entries. Then one cal- 
culates the Jacobians of the change of variables Mi — > (X, ?7i),M2 (X,U2) : dMi = 
{A{ic)fUl=idxkUi<jd^iUi)i,dQ{Ui),, and dM2 = iAiy)f Uk=ldykU^<J d^iU2)^JdQiU2)^J. 
As a result we get 



P^^{n,t,t*) = J eS-iEr=i*'"-reEr=iE:=i*-2'r/(x,y)A2(x)A2(y) f[ dx^dy^ j dMi 

(5.3.3) 



fc=i 



where 



/(x,y)= / e-Tr^^^^-^f/ (5.3.4) 



is the so-called Harish-Chandra-Itzykson-Zuber (HCIZ) integral. The integration ( ^.3.4|) 
is over unitary group U{n), d^U is Haar measure. This integral was calculated : 

det (e-'^'yk^Y 
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Due to ( p.l.24D formula ( p.3.5| ) is a manifest of the fact that the Harish-Chandra-Itzykson- 
Zuber integral actually is a tau-function (see a different approach to this fact in ^^). This is 



the tau-function which corresponds to r(n) = — l/^i, in other words we can write 

^^""'^^ ri/,,(n,t(x),r(y)), l/r'{n) = -- (5.3.6) 



0!---(n-l)! ^"^'^ n 

For n = 1, xi = x, 2/1 = y we have 

rv,,(l,t(x),t*(y)) = e-^^^ (5.3.7) 

This was the subject of the Example (|5.1.4| ). 

Therefore due to (|5.1.26| ) we obtain that the model of two Hermitian random matrices is 
the series of hypergeometric type ( p.2.13|) , where = Y^=\ ^^^1 r{n) = n: 



j^*^^^(n,t,r) = < ei:t.i::^.*-r,ei:L.i::.i*^^r -^EHA^A(t)s.(r) 

(5.3.8) 

where we remember that according to (|2.2.1| ) 



in),- n(r. I 7 .,_ r(n + l + AOr(n + A.)---r(A.) 

^ ^ ,iv ' ^" r(n + i)r(n)...r(i) ^^-^'^^ 

and Cn incorporates the volume of unitary group. 

The same result one can obtain with the help of ( |5.1.21 ), using the known result that integral 
(|5.3.3| ) is equal to 

Cn r ■■■ r eS-iS-=i*'"^^eSr=iE::=i*™?/re-ELi^«2/-^A(x)A(y) J]dxi% (5.3.10) 

J — OO J — OO 

where Xj and j/j are eigenvalues of Mi and M2 respectively. 

This integral has a form of (|5.2.5|) if one takes the integrals over variables Xi be along the real 



axe, while the integrals over variables i/i be along the imaginary axe. Then we take r(n) = n 
so that ^ ^ 

l^rixy) = n/r' = l-xy + ^- + --- = e-"^ (5.3.11) 



e-^Mzdy = 2tt (5.3.12) 
Thus we have the following perturbation series 



2MM(^^ 0,0) EHA^A(t).A(t*) (5.3.13) 



In case all higher coupling constants are equal to (namely, in the case of Gauss matrix integral) 
this series can be easily evaluated as 



) = — 

and for n > 1 one can use the formula (for example see |^) 



00 exp _ * 
h,t2, 0, 0, . . . ; , t;, 0, 0, . . .) = ^ m\hUt)hm{n = " (5.3.14) 



, n-1 

1 / Q^"^""^ 



Tr(n,t,t*) = -det r,. (l,t,t*) (5.3.15) 

' ' ' (n- l)!---0! ydldp^ r v , , ; 1 v ; 



k '112=0 



which is true for the case r(0) = 0. 
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5.4 Hermitian one matrix model 



Below we shall follow 0. 

We suggest a fermionic representation of the well-known partition function of one matrix 
model (which generalizes the known Gaussian unitary ensemble (GUE), ^) 



Z{N,g,g,) = J rfMe-^Tr(fA^2H-34M*)^ 



(5.4.1) 



this fermionic representation being different from pH|,[pD|. The reason is that, as it was before 
in the paper, we consider the matrix model in the sense of the perturbation series . 

We chose the normalization of the matrix integral in such a way that it is equal to 1 when 
94 = 0. 

Let us evaluate simplest perturbation terms for one matrix model using the series (|5.3.14|) 
Let us consider the model of two Hermitian matrices ( |5.3.1| ), and take all tk = except t^, 

3 known |^ if we put 

-AN-%, g = {2Ntl)-^ (5.4.2) 
/2*^^^(0, 0, 0, t4, 0, . . . ; 0, t*, 0,0,...) = Z{N, g, g,) 



and all t^. = except tg- Then as it is known |^ if we put 

94 

we obtain that 



and therefore we get 

Y^{N),s,{t)s,{t*) = J rfMe-^Tr(|M^+,4M^) ^ z{N,g,g,) 



(5.4.3) 



(5.4.4) 



Remark. As it is well-known, see for instance jSQ], according to Feynman rules for one 
matrix model we have (a) to each propagator (double line) is associated a factor l/{Ng) (which 
is in our notations) (b) to each four leg vertex is associated a factor {—Ng4) (which is At^ 
in our notations) (c) to each closed single line is associated a factor A^. Therefore one may say 
that the factors (A^)a in ( ^.4.4|) is responsible for closed lines, the factors sx(t = 0, t2, 0, . . .) are 
responsible for propagators and the factors SA(t*) are responsible for vertices. As we see we get 
the following structure for the Feynman diagrams, containing k = |A|/4 vertices and 2k = |A|/2 
propagators: 

iNg,y'\/\Ngp'\/' ^ {N),axbx (5.4.5) 

|A|=4A: 

where the number |A| is the weight of partition A, and ax, bx are numbers: 



ax = sa(0,0,0,1,0. 



6a = sa(0,1,0,...) 



(5.4.6) 



Let us also notice that to evaluate the numbers ax, bx one can use either the definition of Schur 
functions ( p.l.26|) -( p.l.27|) , or one can use the following formulae, which we borrow in ||56|| : 



7V-1 A(/i'=^'=")A(/i°'^'^ 
11 ^■ 



k=l 
N-1 



h<"id_l 



Sign 



h 



odd\ 



ax 



k=l £=0 Hi ) 



Sign 



n n(/^ 

0<ei<€2<3 i,j 



(5.4.7) 



(5.4.8) 



where we keep the notation of |5g] for shifted weights Xi — i (not miss with the complete 
symmetric function). Here integers h factor into 4 groups of ^ integers h"^ with e = 0,1,2,3 
denoting their congruence modulo 4 ||56|. □ 
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Since only is non vanishing the first non vanishing Schur functions correspond to A of 
weight 4, which are (4), (3, 1), (2, 2), (2, 1, 1), (1, 1, 1, 1). Let us notice that two last partition 
are conjugated to the two first, and the partition (2, 2) is self conjugated. 

To evaluate l.h.s. of ( |5.4.4| ) we use (|2.2.1|) , ( |2.1.12| ) and take into account the property for 



conjugated partition A' we have 

sx'{-t) = (-)l^lsA(t), r'y{-n) = rx{n) (5.4.9) 

where |A| is the number of boxes in the Young diagram (or the same - the weight of the 
partition), and for r' see ( ^.4.91 ). 



The partition A = (4) gives 



N{N + 1){N + 2){N + 3)t 



it 



*\2 
2) 



2! 



{N^ + + IIN^ + 6N) (-^) (5.4.10) 



Due to ( p.4.9|) , ( p.4.2|) the conjugated partition (1, 1, 1, 1) gives the same answer but 



—N. (Thus one needs to keep only even powers of A^). 
The partition (3, 1) gives 

N{N+1){N + 2){N-1){-U)^-^ = {N' + 2N'-N'-2N)[^-^^ (^^^ (5.4.11) 

Again due to the contribution of the conjugated partition (2, 1, 1) one needs to keep only 
even powers of A^. 

First Schur function in the l.h.s of (|5.4.4| ) is equal to zero (due to ( 12.1.121) and since only 



is non vanishing). 

Finely we find first perturbation terms as 



2 ^4] 

This answer coincides with the well-known answer one gets by Feynman diagram method, 
see [|g 

For the next order one takes partitions of weight 8, as a consequence that it is only that 
is non vanishing. One can get that non vanishing contribution is only due to self conjugated 
partitions (3, 3, 2), (4, 2, 1, 1), and to partitions (8); (7, 1), (4, 4); (6, 1, 1), (4, 3, 1); (5, 1, 1, 1) and 
conjugated to these ones. The enumerated partitions yields respectively hi, hi, then hg, —hs, hi; 
—hs, —hi; —hs for ■SA(t*). And it yields respectively hl-h^h^, hl-h^^h^, then h^; -hg, hi; -hg- 
^6^2, —hi + /i6^2; —^8 + ^6^2 for ■SA(t). For partitions (3,3,2), (4,2, 1, 1), and for partitions 
(8); (7, 1), (4, 4); (6, 1, 1), (4, 3, 1); (5, 1, 1, 1) we get respectively 

N{N+l){N + 2)iN-im{N+l)iN-2)iN-l){tl)((^^^ - ^^{tlf] (5.4.13) 
N(N + 1){N + 2)iN + 3)iN- l)iN)iN -2)iN- 3){tl) | ( ^] - ^{tlf 1 (5.4.14) 



and 



2! 



N{N + 1){N + 2){N + 3){N + A){N + 5)iN + 6){N + 7) ( | j ( ^ j (5.4.15) 
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A^(A^ + 1)(A^ + 2)(A^ + 3){N + A){N + 5)(A^ + Q){N - 1) f -M f — %^ ) (5.4.16) 



2! / I 4! 

2 



N{N + 1){N + 2){N + 3){N- 1){N){N + 1){N + 2) (tfj (^^^ (5.4.17) 

N{N + 1){N + 2){N + 3){N + 4)(iV + 5)(iV - 1){N - 2) (^|^ - t^^^ (5.4.18) 

iV(iV + 1){N + 2){N + 3){N- l)(iV)(iV + l)(iV - 2) (-tl) (^^^ + t*2^^ (5.4.19) 

N{N + 1){N + 2){N + 3){N + A){N -1){N- 2){N - 3) (^-|^ + ^2^-) (5-4.20) 
At the highest order at which is A^^ one gets zero. 

(32A^^ + 320A^^ + 488N^)tl{tlY = {32N^ + 320N^ + 488A^2) (^j = (5.4.21) 



%(32A^^ + 320A^2 + 488) 
9 



Finely 



Z{N, g,g,) = l-^(^ + ^\+ M(32iv4 + 320N' + 488) + O (^) (5.4.22) 



5.5 Unitary two-matrix model 

Let us consider the integral 

/^^^(n,t,t*)= / eS:^i*™Trf/reE:.iCTri/,-gTrt/rV-^^f;^^^f;^ 

Ju{n) 

This model was considered in [^. By the Example 3 of the subsection 4.1 we have 

1 [n,t,t)-2_^ ^^^^ -^"A(a)A(6) ^^'^-^^ 



(the last formula we obtain via ( |4.1.11| )).This expression coincides with the series HCIZ integral 



(|5.3.6| ). First this model was identified with HCIZ integral in . 



5.6 New multi-matrix models which can be solved by the method 
of orthogonal polynomials and the Schur function expansion 

Here we consider new solvable multi-matrix integrals which we obtain with the help of the 
subsection Angle integration of tau functions. Integration over complex matrices. 
Let us remind that the multi-matrix model of Hermitian matrices 

J = 1 gyi(M0+-+v^(Af^)gTrMiM,gTrM2Af3 . . . ^Tim^^.m^ ■ ■ ■ dMj, (5.6.1) 
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where Mi, . . . , are Hermitian n by n matrices and 

oo 

Ffc(Mfc) = ^ i(^)l¥Mr, A; = l,...,iV, (5.6.2) 



m=l 



was shown to be multi-component KP tau function and solved by the method of orthogonal 
polynomials. 

Now we consider more general integral over matrices Mi , . . . , M^r 



j e^i(^^i)+-+^^(^^^)i?i(MiM2)i?2(M2M3) ■ ■ ■ i?^_i(Mjv-iMjv)c/Mi ■ ■ ■ (IMn (5.6.3) 

If we choose the interaction term as hypergeometric function of matrix argument (see ( |4.3.1| )) 
i?fc(il/4Mfc+i) = r,Jn,t,MfcM,+i), t = (1, 0, 0, 0, . . .) (5.6.4) 



then due to (|4.3.2D it is possible to perform the angle integration over each Uk{n) where 
Mk = UkX^'^^Uil, X(^) = diag{xf \ . . . ,x^^'>). Really, ( |4.3.2|) is a tau function which has a 



determinant representation (^4.1.11|) . Finely we obtain the following integral over eigenvalues 
xf^ (z = l,...,n;fc = l,...,iV): 



/ = c / \{pk{s^^P)P2{s^:^) ■ ■■pN-i{xf''^4''^) n eS"^iS:^i*-'(^'''')'"rfxf ) (5.6.5) 

1=1 k=l 

where 

Pkix\'^xt'^) = (1, t, xPxt\ t = (1, 0, 0, 0, . . .) (5.6.6) 

(In the case of the familiar multi-matrix model ( |5.6.1|) each pk{xy), k = 1, . . . ,N is e~^'^). 

The integral (|5.6.5| ) may be evaluated by the method of orthogonal polynomials {pm{t^ , x)}, 
{TCmit^ , x)} which depend on the collection of times t^ = (t^^^ t^^^ . . . , t*^^^) 



Pr 



,(t^, x)7in{t^, y)uj{t^, X, y)dxdx = 5nme'^-^'"'\ n, m = 0, 1 , 2, . . . (5.6.7) 



Each of polynomials is of a formpm(t^,x) = Y.n<mPnm{t^)x'^ , 7rm(t^,x) = Y.n<m'^nm{t^)x'^ , 
and the weight function is 

N-l 



uj{f,x,y) = e^^(^)+^-(^) j p,(xx(2))p2(x(2)a;(3)) ■ ■ ■ p^_i{x^''-^'^y) \{ e""^^^'"'^ dx^^^ (5.6.8) 
We obtain 



k=2 



n 



J = c n e'^'^^*^ (5.6.9) 

k=l 

Let us notice that the orthogonal polynomials are Baker- Akhiezer functions of TL hierarchy. 
Gauss-Zakharov-Shabat factorization problem for TL equation directly results from the 
relation (|5.6.7| ) 

ir+(t^) = ir„(t^)e*G(t^), = dtag{(f)i{t'^),Mt^),---) (5.6.10) 

where 

(i^+)„m = (vr"-^)„m, {K-)nm=Pnm, Gnm = i x"^ y'^Uj{t^ , X , y)dxdx (5.6.11) 



(as it was in the case of two- matrix model |}44| ). 

The Schur function expansion for these matrix matrix models is given by ( p.. 0.9 ), where K\^j_ 
are expressed via products of skew Schur functions. 

34 



5.7 Standard scalar product 



Let us suggest a realization of the standard scalar product ( 3.3.2 ). This scalar product is 
a particular case of ( |3.3.9|) . In this case r = 1. This function has no zeroes, therefore the 
consideration of the previous subsection is not valid. We realize ( p.3.2| ) as follows 



< 



/CO 
/(t)e-S::^i-l*-l'<7(t) n 
m=l 



IT 



(5.7.1) 



where is the complex conjugated of t^- The formula for the scalar product of two tau- 
functions is similar to (|5.1.29| ): 



where due to (|3.4.8D,(p.3.2|) r3(i) = ri{k + z)r2(m + i), and 



„ oo 

/ Tr, {k, t, 7) e- Sr^i '=1^'=! (m, 7, t*) J] 

1=1 



kd'^jk 



TT 



(5.7.2) 



(5.7.3) 



Different representation of the standard scalar product one can get in the following way. 
Consider 

<Tr,{k),Tr,{m)>n= f Tr,{k,t,U)Tr,{m,U-\t*)dM = ('"l'^2)ASA(t)sA(t*), (5.7.4) 

JU(n) u\\^ 



l{X)<n 



where the tau functions of matrix arguments are defined by ( [4.3.1 ). In case the restriction 
/(A) < n is irrelevant the standard scalar product coincides with the scalar product <, >„ of 

5.8 About spectral determinants 

In many problems it is interesting to evaluate the following integral over n by tt, matrices 



K{n : a, b, c, d) 
where 



Uti det(a. - Ml) Uti detjh - M,) 
n£i det(Q - M2) nr=i det(rf, - M2) 



(5.8.1) 

00 00 
Vl(Mi) = imTrMf, V2{M2) = ^ CTrM^"^ (5.8.2) 

m=l m=l 

The integral ( ^.8.1|) is called a spectral determinant. In case of one matrix model the results for 



the product of the determinants were obtained in |6^. In the paper |^ the interesting results 
were obtained to evaluate the rational expressions of the determinants. 
In case the integral 



/(t,t* 



MM,)+V,(M,)^UiM„M2)^M,dM2 



(5.8.3) 



is a tau function of hypergeometric type, the spectral determinant is also a tau-function. By 
the Wick theorem it is equal to a determinant. Really K{x,y) = 



(n|e^(*)^(-) . . ■ ■■r{-)9r{bi) ■ ■ ■ ^*(&pj^(rfpj ■ ■ ■ ^(di)e^^*' V) (5-8.4) 



"Pa 



Cl 



Look Appendix A below. The explicit expression is rather long and we shall write down it 
somewhere else. 
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6 Appendices A 

6.1 Baker- Akhiezer functions and bilinear identities [|7|] 

Vertex operators Voo{t,z), V^(t,z) and Vo(t*,z), VQ(t*,z) act on the space C[ti,t2, ■ ■ ■] of 
polynomials in infinitely many variables, and are defined by the formulae: 

VUt,z) = ^A^e«(*'^)e-«(^'^"'), \C(t,z) = ^-Me-«(*'^)e«(^'^"), (6.1.1 

Voir, z) = ;2Me-«(**'^"')e«(^*'^), V*{t*, z) = ^-MgC(t*,.-i)g-?(a*,.)^ ^q ^ 2 

where 8 = (-^ --^ --^ ) d* = (-^ --^ --^ ) 
We have the rules of the bosonization: 

(M + l|e^(*V(^) =Ko(t,^)(M|e^(*), (M - l|e^(*V*(^) = Ki(t, ^)(M|e^(*) — , (6.1.3 

^{z)e''*^'*^\M) = Vo{t\z)e''*^'*^\M -l),r{z)e''*^'*^\M) = V*{t\z)e''*^'*^\M 
The Baker-Akhiezer functions and conjugated Baker-Akhiezer functions are: 

w^{M,tX,z) = y^^^^, <(M,t,r,^) = ^^^^^, (6.1.5 

Voit*,z)T{M + l) ,,,,,,, . Vo*{t*,z)T{M-l) 

Wq{M, t,t ,z) = -p^^ , WoiM, t,t ,z) = -p^^ , (6.1.6 

where 

r(M) = r(M, t,t*) = (M|e^(*)^e^*(**) |M). (6.1.7 
Both KP and TL hierarchies are described by the bilinear identity: 

^ Woo{M, t, t*, z)w*^{M', t', t'*, z)dz = j wo{M, t, t*, z-^)w*o{M' , t', t'*, z-^)z-^dz, (6.1.8 

which holds for any t, t*, t', t'* for any integers M, M'. 

The Schur functions sx{t) are well-known examples of tau-functions which correspond to 
rational solutions of the KP hierarchy. It is known that not any linear combination of Schur 
functions turns to be a KP tau-function, in order to find these combinations one should solve 



bilinear difference equation, see |2^, which is actually a version of discrete Hirota equation. 
Below we shall present KP tau-functions which are infinite series of Schur polynomials, and 
which turn to be known hypergeometric functions ( ^.0.20| ),( |8.0.26| ). 



6.2 Hq{T), fermions iIj{T, z),tlj*{T, z), bosonization rules |]T2 



Throughout this section we assume that r,r ^0. 
Put 

r(n) = e'^"-^"'^", r{n) = e^'^-''^" (6.2.1) 

New variables Tn,n G Z (and also T„) are defined up to a constant independent of n. 
We define a generator Hq(T) G (7/(00) : 

OD 

H,{T) := Tn-. Vn^n :, (6.2.2) 
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which produces the following transformation: 



(6.2.3) 



This is the transformation which sends Am,Am of ( 3. 2.14]) , ( p.5.2 ) to H_m,Hm respectively. 
Therefore 

e^«(^)i(T)e-^°(^) = if(t), e-^°(^M(t*)e^°(^) = -H*it*). (6.2.4) 
It is convenient to consider the following fermionic operators: 



n=+oo 



^(T,z)=e^°(^V(^)e-^"('^)= E e-^"^>n, 

n=— oo 

ra=+oo n=+oo 



(6.2.5) 
(6.2.6) 



Now we consider Hirota-Miwa change of variables: t = ±t(x^) and t* = ±t*(y^). 
Considering (|6X3|) , ( |6X4D and using (|;2j) , (|;2j) , (g]!]!) (|6X|),(|T|) we get the follow- 
ing bosonization rules 



-A(t*(y^)) 



|M) 



-A{-t'{y^)) 



|M) 



^(T,|/i)---^(T,|/^)|M-iV) 
A+{M,N,T,y^) 
r{T,y,)---r{T,yN)\M + N) 



{M\e 

^^|gA(-t(x-)) 

The coefficients are 

A±(M,iV,T,y^) 

A^(M,iV,t,x^) = 



A-(M,iV,T,y^) 
A(T(x-)) _ (M-iV|^-(-t,^)...^-(-t,J^) 



A+(M, A^,T,x^) 



A-(M, A^,T,x 



nfkM^-Vj) ^, r(M, 0,T,0) 



X 



(l/i---2/^)^^*^ t{MtN,0,T,0)' 
nZ:,i^^-XJ) r(M,0,t,0) 



[Xi ■ ■ ■ Xtv 



|Af-l=FM 



r(MTiV,0,T,0) 



(6.2.7) 
(6.2.8) 

(6.2.9) 
(6.2.10) 

(6.2.11) 
(6.2.12) 



where in case = 1 the Vandermond products Hi^jivi — Vj), Hi^Axi — Xj) are replaced by 1. 
The origin of the notation r(M, 0, T, 0) is explained latter, see ( |6 3.5|) , (|6.3.71 ), and used for 



n-l 



r(n,0,T,0) = 



k=0 



r(0,O,T,O) = 1, n = 0, 



-n-l 



\k+n 



n<0. 



(6.2.13) 

(6.2.14) 
(6.2.15) 



r(n,0,T,0) = g-^^--^-^-^- = e"^° [] (K'^)) 

fc=0 

Therefore in Hirota-Miwa variables one can rewrite the vacuum expectation of the exponents 
of ( P^TTD and (PX^ ): 



(M - N\r i-T, ^)... ri-T, ^)V^(T, yi) ■ ■ ■ y;v) |M - N) 



A+(M, iV, T, x^)A+(M, N, T, y^) 
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(6.2.16) 



^^|gA(-t(x-))g-A(-f(y-))|M) = 



A- (M, N, T, x^) A- (M, AT, T, y^) 
This representation is suitable for an application of Wick rule ( |3.1.6D , see ( |6.7.4| ) below 



(6.2.17) 



6.3 Toda lattice tau-function r(n, t, T, t*) [|T2 



Here we consider Toda lattice (TL) |T^ , [0 ' [0 • notations M, t,t* correspond to the nota- 
tions s,x, —y respectively in p7| . 

Let us consider a special type TL tau-function (|3.1.25| ), which depends on the three sets of 
variables t, T, t* and on M & Z: 



r(M,t,T,t*) = (M|e^(*)exp K:T„ : : e^*(**)|M), 



(6.3.1) 



where : ip^i^n '■= i'ni^n — {Oli'ni^nlO) . Since the operator J^'^oo '■ i^ni^n '■ commutes with all 
elements of the gl{oo) algebra, one can put T_i = in (|6.3.1| ). With respect to the KP and the 
TL dynamics the variables T„ have a meaning of integrals of motion. 

As we shall see the hypergeo metric functions ( |8.0.1CI| ),( ^.0.13|) , (|8.0.20|) ,( |8.0.171 ) listed in the 
Appendix are ratios of tau-functions (|6.3.1|) evaluated at special values of times M, t, T, t*. It is 
true only in the case when all parameters of the hypergeometric functions are non integers. 
For the case when at least one of the indices is an integer, we will need a tau-function of an 
open Toda chain which will be considered in the Appendix. 

Tau-function ( 6.3. Ij ) is linear in each e^". It is described by the Proposition 



Proposition 10 

r(M,0,T,0) Av ; Av V • • ; 

The sum is going over all different partitions 

A= (Ai,A2,...,A0, / = 1,2,3,..., (6.3.3) 

excluding the partition 0. 



Let r 7^ 0,r 7^ 0. Then we put 



r(n) = e^"-^ r(n) = e^""^ 



to show the equivalence of ( p.4.8|) and ( |6.3.2| ) in this case. We have 

n-l 

T{n, 0, T, 0) = e-^"-i--^i-^« = e""^" J] (^('^))""'' > ^ > 0. 

fc=0 

r(0,O,T,O) = 1, n = 0. 



-n-l 



Tin, 0, T, 0) = e^"+-+^-2+^-i = e""^" H (K'^))"'"" ^ ^ < 0. 

fc=0 



(6.3.4) 

(6.3.5) 

(6.3.6) 
(6.3.7) 
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Proposition 11 Lei r(n, t, T, t*) is the Toda lattice tau function ( \6. 5'. Then 

r(n,t,T + T,t* 



r(n,0,f + T,0) 



(n|e^(*)e-^(**V) = r,>(n, t,t*). 



where 



A{t*) = J2A,tl A{t) = J2^ktk, 



(6.3.8) 



(6.3.9) 



A;=l 



fc=l 



Ak,Ak are defined by ( ^. 5. !^ ), ^3. 5.4 )- The functions r,r are related to T,T by (\6. 3.4 ), the 
product rf has no zeroes at integer values of it's argument. 
This proposition follows from formulas ( \6. 2.d{ )-( \6.2. 4\ ) ■ 



We shall consider f = 1. 

The notation Tr{M, t,t*) will be used only for the KP tau-function ( p.4.6| ). Notation 
r(M, t,T, t*) is used for Toda lattice tau-function. These tau-functions are related via ( |6.3.8|) . 
KP tau-function T{n) := Xf.(n, t,t*) obeys the following Hirota equation: 



T{n)dqdt^T{n) - dt^T{n)dqr{n) = r{n)T{n - l)r(n + 1). 

The equation 

dtAl(t)n = r{n)e^^-^-'^^ - r{n + l)e<^"-<^"+i 
which is similar to the Toda lattice equation holds for 

7:,(n + l,t,t*) 



/>n(t,t* 



-log- 



r,,(n, t,t*) 



(6.3.10) 
(6.3.11) 

(6.3.12) 



Equations (|6.3.11|) and ( |6.3.10|) are still true in case r has zeroes. 
If the function r has no integer zeroes, using the change of variables 



0n -^n; 



we obtain Toda lattice equation in the standard form p7[] : 



(6.3.13) 



(6.3.14) 



As we see from (|6.3.13|) the variables T„ might have the meaning of asymptotic values of the 
fields 0n for the class of tau-functions ( |6.3.1|) which is characterized by the property as 

ti 0. 



6.4 Expansion of and its trace in the Fock subspace F„ 



n\ 



n=—oo A 



It then follows that the loop 



« r,(n) »r,n= E^A(n) = Tr^„e^«(^) 



(6.4.1) 



(6.4.2) 



where the trace is taken in the n-th component of Fock space Fn 
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6.5 Expressions and linear equations for tau-functions I (Hirota- 
Miwa variables) 

It is well-known fact that tau-functions solves bilinear equations. The class of tau-functions 
under our consideration solves also linear equations. 

The relations of this subsection correspond to the case when at least one set of variables 
(either t or t* or both ones) is substituted with Hirota-Miwa variables. 

In cases t = t(x^) (using (^2^) and (|3A^ )) and t = -t(x^) (using §JJ^) and (|3X^ )) 
one gets the following representations 



r,(M, t(x^), t*) = A-^ (^e«.'(t*.^i) . . . e€^'(**'^^) • a) , A = 
r,(M, -t(x^), t*) = A-^ {e-^r{t',x,) . . . g-c.(t*,x^) . ^\ ^ 



(xi---xjv)^-i-^'^ 

~ N-l+M 



{xi ■■■xn) 



(6.5.1) 



(6.5.2) 



where the exponents are formal series = 1 + ^ + ■ ■ ■, and r{Di) Hk — '"'i.'^i) Ylk ■ Let us 
write down that according to (|3.4.4| ),( p.4.9|) 

+00 +00 Q 

e.Kt*,x,)= ^C(x,r(A))"^, Ui*,x)=J2tl{xir{-D,)r, A = x,— (6.5.3) 

m=l m=l 

Let US note that [C^r'{'t*,Xn),^r'{'t*,Xm)] = = [C,r{'t*,Xn),C,r{'t*,Xm)] foY all n,m. 
( |6.5.2| ) may be also obtained from ( |6.5.1| ) with the help of ( p.4.10| ). 
From (|6.5.1|) it follows for m = 1,2,... that 



d 



N 



dt 



--Y^(x.r{D,jr) Ar,(M,t(x^),r) =0, A 



Jlj^jf (Xj Xj 



rn i=l 



(xi---x^)^-i-A^^' 



In case t = t(x^),t* = t*(y^ ) we get 



r,(M,t(x^),t*(y^*)) = ^r„(M).,(x^)sA(y^*) 

AeP 

1 N N* 



Mx 



nn(i- wpxj)-'-A( 



X 



i=l j=l 

N N* 



YlUil - x,y,r{Dy^))-' ■ A{y) 



where 



A(x) 



Ll<J 



1— r /V* / 



[xi - ■ -xnT-'-'" ■ ■ (yi- ■■?/iv)^*-^-*^ 
Looking at (|6.5.6| ) , (|6.5.7D one derives the following system of linear equations 



(6.5.4) 

(6.5.5) 
(6.5.6) 
(6.5.7) 

(6.5.8) 



- E , , . +n] (A(x)r.(M,t(x^),t-(y^*))) = 0, j = l,...,iV* (6.5.9) 



N* 



1 



iV*|(A(y)r,(M,t(x^),r(y^^)))=0, ^=l,...,iV (6.5.10) 



In some examples below we shall take specify t* as follows 
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(1) Let us take 



N 



mtm = Y.^k^ r = (1,0,0, 

k=l 



Then one consider the equation (|6.5.19D with k = 1 (taking t* = 1) 



where 



and Di = Xid^^- We get 



r,(M,t(x^),t*) = _e^iriD,)+-+x^r{D^) . 



A 



E A - ^ (j^x^riD,)^ Aj r,(M,t(x^),r) = 



(6.5.11) 



(6.5.12) 



(6.5.13) 



(6.5.14) 



(2) Let us take 



N 



fc=l y 



m = 1,2, ... 



(6.5.15) 



Let us note that the tau-function rr(t,t*) will not change if instead of ( |6.5.15|) one take 

1 



N 



mtm = J^i^^k)""^ "^C = 1 

k=l 1 - ^ 



■, m = l,2,... 



(6.5.16) 



This gives the representation as follows 

r,(M,t(x^),r) = 

We see that 



N 

n 



A(x) {xir{D^J;q)^ 



A(x) 



(6.5.17) 



(6.5.18) 



i=l 



We shall write down more linear equations, which follow from the explicit fermionic repre- 
sentation of the tau-function ( |3.4.6|) via the bosonization formulae ( 6.2. 16| ) and ( 6.2. 17| ) These 
equations may be also viewed as the constraint which result in the string equations (|6.14.4|) 
and ( |m:5|) . 

For the variables in case t = t(x^), one can use the relation {M\Am = 0, and makes 
profit of the relation = e^°H_me~^° inside the fermionic vacuum expectation value ( [6.2.161 ). 
This way we get the partial differential equations for the tau-function (|3.4.8|) : 



ar,(M,t(x^),r^ 
dt* 



i (j^{x.,r{D,,;)r^ ATr{M, t(x^), r ), m = 1, 2, 3, . . . , (6.5.19) 



where A is proportional to A+ of ( 6.2.12 ): 



A 



(xi---x^)^-i-A'^' 



(6.5.20) 
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In variables y*-""-*, t* = — t*(y(°°)), we can rewrite ( |6.5.19| ): 



E 

i=l 



dyi 



N 



(6.5.21) 



where efc(y) is a symmetric function defined through the relation nS(l + %) = Z^fciS ^^^^(y). 
Also we have 

/M+N-l N \ _ 

E fc-E^:^0^"^(A^'-*^^'')'^'-**(y^'''^))^" = 
\ fc=i i=i / 



' M+N'-l 



N' 



E ^ - E (^-A.,?/. ) 1 A-r(M, -t(x^),T, -r(y(^')))A- 



(6.5.22) 



fc=i 



where A" = A-(M, A^, 0, x^) and A" = A-(M, A^, 0, y(^')). This formula is obtained by 
the insertion of the fermionic operator res^ '■ ip*{z)z-^ip{z) : inside the fermionic vacuum 
expectation. 

These formulae can be also written in terms of higher KP and TL times, with the help of 
vertex operator action, see the Subsection "Vertex operator action". Then the relations (|6.5.19|) 
are the infinitesimal version of ( |6.6.7 ), while the relation (|6.5.22| ) is the infinitesimal version of 
(ip). 

In some examples below we shall take specify t* as follows. (1) Let us take 



mtr. 



N 

E 

fc=i 



X 



k 5 



(tt,o,o,...) 



(6.5.23) 



One notes that that the tau-function Tr{t,t*) will not change if instead of (|6.5.15|) one take 

TV 

mt^ = Y.{tiXkT, t* = (1,0,0,...) (6.5.24) 

k=l 

Then one consider the equation (|6.5.19| ) with = 1, taking ti = 1. We put Di = Xid^^ and get 

(E A - ^ (E A)) r,(t(x^), r ) = (6.5.25) 
where r is given by (|6.12.1|) , and 

ni<j {Xi Xj^ 



A 



[xi--- Xn 



\N-M-1 



(6.5.26) 



6.6 The vertex operator action. Linear equations for tau-functions 
II 

Now we present relations between hypergeometric functions which follow from the soliton the- 
ory, for instance see |T^,[^. (In there were few misprints in the formulae presented below). 
Let us introduce the operators which act on functions of t variables: 



ni^\t, t*) := - limres,=oC(t, ^ + z-')V^{t, z 



dz 



(6.6.1) 



42 



dz 

^]W(t^t) :=limres,=oV^o*(t*,z + e)e(°Ht,^)K)(t*,^)-, 
where Voo(t, z), V^(t, z), Vo(t*, z), lo*(t*, z) are defined by ( p.l.l| ), and 

• -j^ ^ m +00 



For instance 



Z] 



m=l ^ ^ 'm=l 



r = 1 : a(°^)(t,t*) = fil°)(t*,t) = E^^nC- 

n>0 

Also we consider 

dzdz* 



Znn{t) = res2=ores^*=o— K^(t, z*)Voc{t, z) 



zz"* 



Z*"' zz* 

The bosonization formulae (|6.1.3| ),( ^.1.4|) result in 



Proposition 12 We have shift argument formulae for the tau function l \3.4-^ : 

gn(°°'(t,7) . ^^^^^ t*) = Tr{M, t, t* + 7), e^'-''(**'^) ■ Tr{M, t, t*) = r,(M, t + 7, t*) 
Also we have 

eE-^7n + M^nn(t) _ ^^^^ ^ rj. ^ ^ ^^^^^ 7n + M (t * ) _ ^^^^ ^ rj. ^ ^ ^ ^ t, T + 7, t* 

For instance 
Also 



6.6.2) 

6.6.3) 
6.6.4) 

6.6.5) 
6.6.6) 



6.6.7) 

6.6.8) 
6.6.9) 



e^-- ^"+*^^"" • exp E ^''^ntn = e^-^ ^"+«^™ ■ exp ^ ^^"^n = ^(M, t, T, t*). (6.6.10) 



\n=l / \n=l 



6.7 Determinant formulae I (Hirota-Miwa variables) [[T2|] 

In the case t = t(x^) one can apply Wick theorem to obtain a determinant formulae. 
Proposition 13 A generalization of Milne's determinant formula 

det (xf-^JM - k + l,t(xi),t*))^ 

det fxf-'^)'^ 

Proof follows from the representation ( p. 5.1 



r,(M,t(x^),t*) = . ..^ ' ^' (6.7.1) 



r,(M,t(x^),t*) = (e«^'{t*--i) . ..e^A^*.-N) . a) , A = ^""f \, , (6.7.2) 

from the explicit form of the Vandermond determinant: Yli^j{xi — Xj) = det(xf'~^) and from 
the representation ( |6.7.2| ) for the case when the number of variables x is equal to one as: 

Tr{M -k + 1, t{xi),t*) = e«'-'(**'^>) ■ xf^-^+^ (6.7.3) 



2 3 

where t(xj) = (x^, ^, ^, . . .). 



In the case t = t(x^),t* = t*(y^) we apply to (|6.2.16 ) the Wick theorem to obtain a 
different Proposition 
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Proposition 14 For r ^ we take a tau function r,.(M, t(x^), t*(y^)) and apply Wick's 
theorem. We get the determinant formula: 

ca4 det (rr(M — + 1, Xj, w,))^. -, 

r.(M.t(x").f(y-))^-ii- \(^.)^(^.;'''^'-' . (6.7.4) 

where 

n-l 

cn = n (KA:))'"" (6.7.5) 

fc=0 

Tr{M-N + 1, Xi, Vj) = 1 + r{M -N + l)xiyj + r{M -N + l)r(M - + 2)x^?/| + ■ ■ ■ (6.7.6) 
One can prove it also either from the representation ( |6.5.6|) or from the representation 



6.8 Determinant formulae II (The case r(0) = 0.) 

There is also determinant formulae for tau functions in variables t, t*. These one we get in case 
the function r has zeroes. If r(0) =0 and n > then 

/ ga+b 

Tr{n, t, r ) = On det — — ^r,(l, t, t* ) (6.8.1) 

\CtlCh J a,b=0 

Cn = n (K^))'^" (6.8.2) 

fc=0 

6.9 Integral representations I (Hirota-Miwa variables) 

In the case t = t(x^), t* = t*(y^) we get an integral representation formulae. For the fermions 



(|6.2.5| ) we easily get the relations: 



Jij{T,az)dfi{a) = ij{T + T{fi),z), J ^j* (^-T, -^^ dfi{a) = ^j* (^-T - T{fi),^^ (6.9.1) 

where /x, /i are some integration measures, and shifts of times T„ are defined in terms of the 
moments: 

J a"d/i(a) = e-^"^''\ J a^'dfi^a) = e-^"^^^\ (6.9.2) 
Therefore thanks to the bosonization formulae (|6.2.16| ) we have the relations for the tau- 



function; below t* is defined via Hirota-Miwa variables. 
Proposition 15 Integral representation formula holds 



r r(M,t(ax^),T + T,t(ay(^))) tv 
/ At(ax ) n ^ , ^ Ariay^^^) [[dfi{ai) [[dfi{ai 



(Af , 0, T + f , 



r M,t(x^),T + T + T(A)+T(;.),t(y^)^ ^ 



r (M,0,T + T + T(/i) +T(/i),0 
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where AT(«y(;v)) = A+(M, iV, T, ay^j.^), A^(«x^) = A+(M, N, f , «x^), 

'^y(Af) ~ (o^il/i; ^21/2, • • • , Q-nVn) o,nd oix.^ = {aiXi, 02X2, ■ ■ ■ , Q-nXn)- In particular 

.r(M,t,T,t(ay(^))) , , 

I riM,0,T,0) ^Tiay,^))nMa.) = 



r(M,t,T + T(/i),t(y^ 
r(M,0,T + T(/i),0) 



AT+TM(y'')- (6.9.4) 



Remember that arbitrary linear combination of tau-functions is not a tau-function. Formulae 
(|6.9.3| ) and also (|6.9.4| ) give the integral representations for the tau-function ( |3.4.6|) . It may 
help to express a tau-function with the help of a more simple one. 
We shall use the following integrals 



27r 



/ tP*(-T,-)e''a-^da = tP*(-T ~T\-), (6.9.5) 

a/— 1 Jc X X 

roo 1 1 

/ ^/^*(-T,— )e-"a'^-Ma = V*(-T-T",-), (6.9.6) 
Jo ax X 

— ^ r ^*(_T, — )a-i(l - a)^— Ma = ^*(-T - 1), (6.9.7) 
L [0 — a) Jo ax X 

where C in ( |6.9.5| ) starts at —00 on the real axis, circles the origin in the counterclockwise 
direction and returns to the starting point, and 

T„^ = lnr(6 + n + l), = - lnr(a + n + 1), = In ^^^tnr//.^ v (6-9-8) 

1 (a + n + 1)1 (0 — a) 

Let us remind that q- versions of exponential function, of gamma function and of beta function 



are known as follows 25 



exp,(x) = , \ , , r,(x) = (l-g)i-/^, i?,(x,^) = y;^yf (6.9.9) 
(a;(l - g);g)oo (g"^; g)oo Tq{x + y) 

Gamma and beta functions have the g-integral representations 

r,(x) = -g-^ r -^^^d,a, Bqix, y) = f a^-'j^^l^d,a (6.9.10) 
Jo exp^(a) Jo {aqy;q)oo 



where g-integral are defined as 

j=oo ]=oo 

E ^W)' / f{x)d,x = c{l-q) E 
The following g-integrals will be of use: 



/ /(x)ci,x = (l-g) E ^W)' f{x)d,x = c{l-q)J2q^f{cq') (6.9.11) 



^ " r ^*(^T,—^-^)^^d,a = r (-T -T{a,q),-) , (6.9.12) 
1 JO a(l — q)x expg(Q;) V xj 



q — \ Jo ' a{l — q)x exp (a 



^ l^—. f'r{-T,—)a--'j^^^^P^d,a = r{-T-T{a,b,q),-). (6.9.13) 
1 q[o — a) Jo ax {aq° ; gjoo a; 
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where, as one can check using ( |6.9.10D , 



Tn(a,q)= In- — 4 r, TJa,b,q) = In + ^ + (6.9.14) 

^ '^^ (l-g)"r,(a + n + l)' , Yg{a + n + l) ^ ' 



In the same way one can consider Hirota-Miwa variables (|1.0.8| ). In the Examples below we 
shall present hypergeometric functions listed in the Subsections 1.2 and 1.3 a.s tau-functions of 
the type ( |6.2.16|) . Then we are able to write down integration formulae, namely (|6.9.4| ), which 



express p+i$s and p+i$s+i in terms of with the help of ( |6.9.12|) , ( |6.9.13|) and ( |6.9.3|) . In 



38 1 different integral representation formula was presented, which was based on the g-analog 



of Selberg's integral of Askey and Kadell. By taking the limit q ^ 1 one can consider functions 
pFg. Using ( |6.9.8| ), one can express p+iJ-'s, p+iJ^s+i and pJ-'s+i as integrals of pjF, with the help 
of (HH), ( |6X7|) and ( |6X5D respectively. 

6.10 Integral representations II 

Now let us consider the case when depends on T,t*. Using (|3.4.7| ) and ( ^.1.21|) one gets the 
integral representation for in case the function r has zero. Let us again choose r(0) = 0. The 
have the representation described above 

Proposition 16 

r,(M,t,r) =/■■■/ A(^)A(^*) n e^'^-A^"^'-^^T^*-) iir{zkz\)dzj,d4 (6.10.1) 
•' k=i 

where 

M M 

H^) = m^^-z,), A(z*)=n(4-4) (6-10.2) 

i<j i<j 



an 



d Hr is defined by l \5.1.1(\ ),( \5.1.3i ). 



Now it is interesting to notice that if the series Ti/r (which defines fir, see ( ^.1.10|) ) is a 
convergent series then the series Tr may be divergent one. In this case one can consider the 
integral /^(M, t,t*) is an analog of Borel sum for the divergent series r^, see 0. 

6.11 Rational r and their deformations 

Lemma 1 If 

. ^ nLl(^ + «i) 11 IX 

r[n) = — — -. — (6.11.1) 



then 



where 



taH = 7 +00 1 , x\ 6.11.2 

nfc=i SA(7(0fe + ^)) 

7(00) = (1,0,0,...), 7(a)= (6.11.3) 



and sa(7) are Schur functions as functions of variables 7 according to (|2.1.9D and (|2.1.12| ).D 
Lemma 2 If 

nLi(i + g"+"') ..11 

r[n) = — - — -TT 6.11.4 
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then 



where 



rxin) = (sA(7(+oo,g))) 



1 - (q") 

7m(a, q) = — li+oo, q) 



m(l — g™) 

□ For the proof we use the following relations (see 0): 



s-p nLi sx{n{ak + n,q)) 
riLi sx{-i(hk + n,q)) 

— -, m = l,2, 

m(l — q^) 



(6.11.5) 
(6.11.6) 



(M)GA SA(7(+oo,g)) sa(7(+oo) 

Now we see that for arbitrary complex h we have 
Proposition 17 

< e<^ , eS™=i >,.,n=< 6^''+") S"=i , eS™=i >^,,„, rf,(fc) : 



r(A;) 
b + k' 



Also we have 



< eZ^m=l l-,™ e^™=l >rr7 = < 6^™= 



°° 7m(l-g' 



1-9" 



(6.11.8) 



r(A;) 

1 - qb+^ ■ 
(6.11.9) 



Let us note that parameters 7m (o, g) and 7m (o) are chosen via 'generalized Hirota-Miwa trans- 
forms' ( |6.11.6| ),( |6.11.3| ) with 'multiplicity a'. Remember that |g| < 1. Due to 



Sx{l{+oo,q)) = lim SA(7(a,g)) 



q 



n(A) 



sa(7(+oo)) = lim sx 

a — >+oo 



'71(a) 72(a) 



Hx{qy 
1 



lim ^SA(7(a)) „ 



(6.11.10) 
(6.11.11) 



This allows to rewrite the series ( |6.12.7|) and ( |6.12.25|) only in terms of Schur functions, see 



12 



For functions pF, of ( |6.12.7| ),( |6J2:9| ) and of ( |6.12.25| ),( |02:26| ) we use the results of 
subsection 5.4 to derive the following representations 



p-^ s 



'ai + M, . . . ,ap + M 
61 + M, . . . , 6, + M 



A(x) 

N 



^x\r{D-i)-\ \-XMr{DM) 



A(x) 



(6.11.12) 
(6.11.13) 



/ai + M, . . . , Op + M 
^ M 61 + M, . . . , 6, + M 



N 

n 



A(x 



A(x) (x»r(A,);g)c 

1 (a;igV(A)(l-g'+''')"';9 
11 



i=l 



{xir{Di);q)c 



■A{x) 



A(x) 



(6.11.14) 
(6.11.15) 



where Di = Xi-^, A(x) is defined in ( |6.5.8| ), the notation (fe; g)oo see in ( |8.0.14| ). 
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6.12 Examples of hyper geometric series 



The main point of this subsection is the observation that if r{D) is a rational function of D 
then Tr is a hypergeometric series. If r{D) is a rational function of we obtain g-deformed 
hypergeometric series. Concerning the hypergeometric functions see Appendix C. Now let us 
consider various r{D). 



Example 1 Let all parameters bk be non integers. 

{D + ai){D + a2) ■ ■ ■ (D + ttp) 



oTsiD) 



iD + b,){D + b2)---iD + bs) 
For the vacuum expectation value ( |3.4.8D we have: 

^rAM,t,t ) = E^A(t)..(t )(^^ + M),...(6. + M),- 
Let us notice that if we put 

tl = l, t* = 0, i> 1, 
then SA(t*) = H^^ (see ( |6.11.11|) ), where Hx is the following hook product: 

Hx = II ha, hij = {Xi + X'j - i - j + 1). 



(6.12.1) 



!A = n ^ii' 



(6.12.2) 



(6.12.3) 



(6.12.4) 



Thus for t* as in ( |6.12.3| ) the tau function ( |6.12.2[ ) yields the hypergeometric function related 
to Schur functions as follows: 



Pr;(M, t,t* 



s 



til ^2 



, >^Z; • • • 



fli + M, . . . , Op + M 
bi + M,...,bs + M 

(ai + M)a • ■ ■ (flp + M)x sxit) 
^ ib^ + M)x--- {bs + M)x Hx ■ 



(6.12.5) 



As an example taking ai = 0, t = (ti, 0, 0, . . .), t* = {t\, 0, 0, 
ordinary hypergeometric function of one variable: 



we obtain the well-known 



p-i 



F,(a2 ± 1, . . . , s ± 1; 6i ± 1, . . . , 6, ± 1; i^*) = r,,(±l, t, t*), r{n) 



n 



m=i{h + n) 

(6.12.6) 



Now we take t = t(x^),t* = (1,0,0,...) then the ( |6.12.5| ) turns out to be the known 
hypergeometric function of matrix argument X 



X 



E 

A 

1{\)<N 



fai + M,...,ap + M 
" '\bi + M,...X + M , 

(ai + M)A---(ap + M)A5A(x^) 
(6i + M)a---(&. + M), Ex ' 



v;(M,t(x^),r) 



(6.12.7) 



N 



This function depends only on the eigenvalues of the matrix X which are x 
The hypergeometric function of matrix argument is also known as the hypergeometric function 
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(|8.0.26|) related to zonal polynomials for the symmetric space GL{N,C)/U{N) 



|. Here 

1, are the eigenvalues of the matrix X, and for zonal spherical polynomials 



there is the following matrix integral representation 



Za(X) = ZxilN) I {U*X.U) dM. 

JU(N,C) 



(6.12.8) 



■ -A^^ and Ai, . . . A^v are main minors of the matrix X, d^:U 



where A^ (X) = Ai^^^^Aa'"^^ 

is the invariant measure on U{N,C), see In our case Z\(X.) coincides with the Schur 

5a 



function s\(:s.^) 



Due to ( |6.5.1|) hypergeometric function ( 5.12.5 ) has the representation as follows 



s 



where Di 



^ai + M,. . . ,ap + M 
^bi + M,...,bs + M 

Xi and 



A 



— exp(a;ir(Di)) • ■ ■ exp(x7vr(L'iv)) ■ A, 



(6.12.9) 



(6.12.10) 



(See also ( |6.11.13| ) for different representations). 

Now let us take t* = (t^, 0, 0, . . .). Then we get the same function (|6.12.7|) , but with each 
Xi changed by —t\xi. Then one consider the equation ( |6.5.19| ) with k = 1. We put Di = Xid^^ 
and get 



fli + M, . . . , Op + M 
+ M,...,bs + M 







(6.12.11) 



where r is given by (|6.12.1|) , and A see in ( 6.12.101 ). 

Taking N = l,t* = (1,0,0,...) we obtain the ordinary hypergeometric function of one 



variable x = Xi, which is (compare with (|6.12.6| )): 
pFs {ai + M, . . . ,ap + M;bi + M, . . . ,bs + M; x) 



-M xr(D) M 

X e ^ ' ■ X , 



D = x4- (6.12.12) 

dx 



The ordinary hypergeometric series satisfies the known hypergeometric equation (we put M 
0) written in the following form: 



(dx - pTsiD)) pFs{ai, . . . ,ap;bi, . . . ,bs;x) = 0, D := xd^. 



(6.12.13) 



This relation helps us to understand a meaning of the function r. 

It is known that the series ( |6.12.7| ) diverges if p > s + 1 (until any of + M is non positive 
integer). In case p = s + 1 it converges in certain domain in the vicinity of = 0. For p < s + 1 
the series ( |6.12.7]) converges for all x^. These known facts (see ||25[) can be also obtained with 
the help of the determinant representation ( |6.7.1|) and properties of (|8.0.10|) . 



Now we have the following representations: 

p+iFs 



{ai ■ ■ ■ e 



-ON TP 
p^ s 





■ ■ , CLp 




h,.. 


■,bs 




Xi 


Oi, . . 


, Op 






bi,.. 


,bs 


aiXi 



a; AT 



U^<j{Xi-Xj) Jo 

N 



J W^Q-iXi — ajXj)dai 
) i<j 



(6.12.14) 

dajsi 
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Hankel type of representation 

ai 



, . . . , Clp 



b,bi, ... ,6, 



c 



(6.12.15) 



where C starts at — oo on the real axis, circles the origin in the counterclockwise direction and 
returns to the starting point. 



Beta type representation: 



/ a,ai, . 


■ 1 (^p 






\bM... 


■ ,bs 







N 



s 



ai, . . . , Op 
bi,...,bs 



aiXi, . . . , aNXN 



N 

\ N N 

] YiiaiXi - ajXj) Yl af^i 

i<j 1=1 



1 — a 



1 

JO 

b-a-l 



(6.12.16) 



Let us notice that according to [|I4| the logarithmic derivative of the hypergeometric function 
(|6.12.5| ) u{x) = log2-Fi (j^ ti = x, 0, 0, . . is a solution of the Painleve V equation. 

Example 2. Let us consider Nhj N matrices X, Y with eigenvalues x^ = (xi, . . . , xat), = 
{Vi, ■ . ■ ^Un) respectively. In order to get a hypergeometric function of two matrix variables 
X, Y, see (|8.0.26|) in the Appendix C, we put 



nr=i(ai + n) 



I{UiiPi + n)N -M + n 



(6.12.17) 



Notice that r now explicitly depends on N . 

Taking into account the relations ( |6.11.7| ) and (|6.11.11|) we get 



(iV)A = //ASA(t(l^)), 1^=(1,1,...,1) (6.12.18) 
and putting t = t(x^) and t* = t*(y^) we obtain the hypergeometric function of ( |8.0.26| ): 



E 

A 

1(\)<N 



X,Y 



ai + M, . . . , Op + M 
6i + M, . . . , 6, + M 

gA(x^)gA(y^) (ai + M)a ■ ■ ■ {ap + M)a 
(iV), (6, + m)a- ■■(&. + M),- 



(6.12.19) 



Due to (|6.5.8| ) hypergeometric function ( |6.12.19| ) has the representation as follows (r = ^r^): 



/ ai + M, . . . , ap + M 
^ M 6i + M, . . . , 6, + M 



N 



(6.12.20) 



Here and below the inverse operator (1 — yjXir{Dx-)) ^ is a formal series 1 + yjXir{Dx-) + ■ ■ ■. 

For example when iV = 1 we get the ordinary hypergeometric function of one variable 
xy = Xiyi: 



,Fs{ai + M,...,ap + M;bi + M,...,bs + M; xy) = x'^' (1 - xyr{D^)y' ■ x*^ (6.12.21) 
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Looking at 12.91 ) one derives the following linear equation 

/ ^ 



D 



i=l 



yjXir{D^ 



+ N 



ai + M,...,ap + M 
bi + M,...,b, + M 



(6.12.22) 



Example 3 The g-generalization of the Example 3 is Milne's hypergeometric function of 
the single set of variables. To get it we choose 

nLi(i-g° 



,rW(n) 



We take t = t(x^),t* = t*(y°°). By putting 



Vk = q 



.k-l U _ 1 f* 



+00 m 



1 



m = 1, 2, . . . 



(6.12.23) 



(6.12.24) 



m m(l — q"^j 

we have SA(t*) = q~'^'^^^H\{q), see ( |6.11.10| ). We get Milne's hypergeometric function (|8. 0.201 ): 



(M| 



^H{t)-A{t*) 



\M) = 



'ai + M, ...,ap + M 
6i + M, . . . , 6, + M 



^ (g^i+*^;g),...(gf^+M;g),i7,(g)''^'' ^• 



(6.12.25) 



A 

i(A)<]V 



According to ( |6.5.8| ),( |6.12.24| ) we have the following representation 



'ai + M,...,ap + M 
bi + M,...,b, + M 



N 

n 



■ A(x) 



= —expg{xir{Di)) ■ ■ ■ expg{xNr{DN)) ■ A 



(6.12.26) 



(the notations (6; g)oo, expg(a) see in ( p.0.14| ), ( |6.9.9| ), A is the same as in ( |6.12.10| )). (See also 
(|6.11.15|) for different representation). 

For this hypergeometric functions we have the linear equation ( |6.5.18| ) (where we choose r 
due to ( |6. 12.231) ), which is g-difference equation. For instance for the case iV = 1,M = we 
get (compare it with ( |6.12.13| )) 

^ (l - g^) - prf\D)^ p$,(ai, . . . , a^; fei, . . . , 6,; g, x) = 0, D := xd„ (6.12.27) 

Milne-s function itself for iV = 1 is the ordinary basic 



where pr'f\D) is defined by 
hypergeometric function: 

' ai + M,. . . ,ap + M 
bi + M,...,bs + M 



q,x 



^ (nbl 



g)n---(g"''+*';g) 



X 



A. (^f>,+A/.^)^...(^f,.+M.g)^ (g;g), 



X = xi(6.12.28) 



Let us take M = 0. We have the following integral representation 

1 



-1$ 



, . . . , U,p 



a, ai 

, bi,...,bs 



q,xi, 



,Xn ] = 



N 



(1-g)^ nf<,(x.-x,)Jo 
q;i(1 — g)xi aAr(l — g)xAr 



(6.12.29) 



n 



dqOtl 
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Also we have the Beta-type representation: 



Qi J Qi'^ ^ • • • ^ 

b,bi,...,bs 



N 



1 /-i 



N 



i<j 



Oil, ■ ■ ■ 1 CLp 

bi,...,b. 



\Ty{b-a)J Ui<j{xi - Xj) Jo Jo 



(6.12.30) 



i=l 



yb-a. n\ ' 

/oo 



Example 4. To obtain Milne's hypergeometric function of two sets of variables x^, y^, we 
take t = t(x^) and t* = t*(y^). This choice restricts the sum over partitions A with Z(A) < N. 
We put 



1 nLi(l-9)"r,(a, + n+l) 

1 - q)^T,{n + N-M + 1) n-=i(l - 9)^5(6^ + n + l)' 



r,(a) = (1 - g) 



Here rg(a) is a g-deformed Gamma-function 



T,{a) = {l-qy 



(g^g)„ = (l-g) 



rg(a) 



„r,(a + n) 



(6.12.31) 

(6.12.32) 
(6.12.33) 



Tq{a) 



Using (PTTTTI ) and ( |6.11.l0[ ) (see III of §) we get 

(g^-*';g)A = g-"('^i/A(g)^A(i,g,...,g^-^) 



(6.12.34) 



(6.12.35) 



(see Appendix C for the notations) and we obtain the Milne's formula ( ^.0.27|) 

' ai + M,. . . ,ap + M 



r.(M,t(x^),r(y^))=,$ 
E 



g,x^,y^ 



61 + M, . . . , 6, + M 



A 

i(A)<iV 



i7A(g) sa(1, g, . . . , g^-i) (g^^^^^; g)A ■ ■ ■ (g^=+*^; q)x 



(6.12.36) 



This is the KP tau-function (but not the TL one because ( |6.12.31| ) depends on TL variable M). 

We have the same type of representation (|6. 12.201) and the same form of linear equation 
(16.12.221) , however we shall use (|6.12.31|) for r in formulae ( |6.12.20|) , (|6.12.22|) . 

To receive the basic hypergeometric function of one set of variables we must put indetermi- 
nates in ( |6.12.25|) as yi = i = {1, . . . ,N). Thus we have 



p^s 



'ai + M,. . . ,ap + M 
61 + M, . . . , 6, + M 



g,x 



N 



(g"i+A^;g);,...(g'^p+M;g);^ g-W 

i(A)<JV 

And for = 1 we have the ordinary g-deformed hypergeometric function: 



(6.12.37) 



'ai + M,...,ap + M 
bi + M,...,bs + M 



\ +g; (g"^+^^;g).---(g"^+"'^;g)n (xy)" rfiioQ«^ 
q,x,y\=^ ,^,^^^f..^ r .^ > ""V = ^1^1 (6-12.38) 



n=0 



{q'^+^';q)n---{q''^+^';q)n {q;q)n 
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6.13 Baker- Akhiezer functions and the elements of Sato Grassman- 
nian related to the tau function of hypergeometric type 

This subsection is based on [12]. Let us write down the expression for Baker- Akhiezer functions 
(|6.1.5| ) in terms of Hirota-Miwa variables ( |1.0.8|) : 

w^iM, -t(x^),r, -) = Ylil - x^+i = (xi,...,x;v,^), (6.13.1) 



r^(M,-t(x^),t* 



i=l 



1. r,(M,t(x^) + [z],t 



N 



We see that the variables Xk,k = 1, . . . ,N are zeroes of Woc{z) and poles of wl^{z). 



dz 



z 



(6.13.2) 



Remark 1 The associated linear problems for Baker- Akhiezer functions are read as 

{dt^ - dt^(f)n) w{n, t, t*, z) = w{n + 1, t, t*, z), (6.13.3) 

dtiw{n, t, r, z) = r{n)e^"~'-^"w{n - 1, t, t*, z). (6.13.4) 

where w is either w^o or Wq. The compatibility of these equations gives rise to the equation 
i\6. 3.11\ ). Taking into account the second eq. ( \6. 3. 1^ ), equations \6.13.^ ), l \6.13.4 ) 'may be also 
viewed as the recurrent equations for the tau-functions which depend on different number of 
variables x^. 

Let us write down a plane of Baker- Akhiezer functions (|6.1.5|) , which characterizes elements 
of Sato Grassmannian related to the tau-function ( |6.3.1| ) rr(M, t,t*). We take = 0, = 
1, . . . , iV in (|6l3l| ) and obt am: 



Woo(n, 0, r, z) = ^"(1+ ^ r(n)r(n-l) ■ ■ ■ r(n-m+l)/i„(-t*)2-'"), n = M, M+1, M+2, .... 

m=l 

(6.13.5) 

The dual plane is 



w*^{n, 0, r, z) = z-'^{l+ r{n)r{n+l) ■ ■ ■ r{n+m-l)h^{t*)z-'^)dz, n = M, M+1, M+2, .... 

m=l 

(6.13.6) 

About these formulae see also ( |6.14.22| ),( |H.14.23| ). 

We see that when r has zeroes, then the element of the Grassmannian defining the solution 
corresponds to the finite-dimensional Grassmannian. 
If r(0) = one can construct a plane 



{<(n,0,t*,l/z),n = l,...,M} 



(6.13.7) 



which is a spanned by a finite number of vectors. The space spanned by these vectors and in 
addition by the vectors ( p. 13.51) gives the half infinite plane generated by {z^, n > 0}. 



53 



6.14 Gauss factorization problem, additional symmetries, string equa- 
tions and ^DO on the circle [[121 



1. We shall also 



Let us describe relevant string equations following Takasaki and Takebe [31 
consider this topic in a more detailed paper. 

Let us introduce infinite matrices to describe KP and TL fiows and symmetries, see ^ . 
Zakharov-Shabat dressing matrices are K and K. K is a lower triangular matrix with unit 
main diagonal: {K)ii = 1. A' is an upper triangular matrix. The matrices K,K depend on 
parameters M, t,T,t*. The matrices (A)jfc = 5j,fc-i, (A)ifc = 5i^k+i- For each value of t,T, t* 
and M & Z they solve Gauss (Riemann-Hilbert) factorization problem for infinite matrices: 

K = KG{M, t, T, t*), G{M, t, T, t*) = exp (^(t, A)) A^^G(0, T, 0)A^ exp (^(t*, A)) .(6.14.1) 



We put \og{Kii) = (f)i+M, and a set of fields 0j(t,t*), (— oo < i < 
higher two-dimensional TL equations. 

^ Take L = KAK~\ L = KAK-\ and (A),^ = idi^k, M = 

M = KAR-^ + M + Y.ntlU\ Then the KP additional symmetries |3l|] , [0 , Q , [|l 
higher TL fiows are written as 



oo) solves the hierarchy of 

KAK-^ + M + EntnL"", 
5i and 



((r(M)L-7" 



d,,K 



dt*K = 

''n 

dt*K 



((r(M)L-i 
-{L-)K. 



(6.14.2) 
(6.14.3) 



Then the string equations are 

LL = r{M), (6.14.4) 

M = M. (6.14.5) 

The first equation is a manifestation of the fact that the group time of the additional 
symmetry of KP can be identified with the Toda lattice time t^. In terms of tau-function we 
have the equation (|6.6.7| ) in terms of vertex operator action |Q,[|l^, or the equations (|6.5.19|) 
in case the tau-function is written in Hirota-Miwa variables. 

The second string equation ( |6.14.5| ) is related to the symmetry of our tau-functions with respect 
to t <^ t*. 
When 

r{M) = M + a, (6.14.6) 

the equations ( |6.14.4| ),( p.l4.5[ ) describe c = 1 string , see [^,|Q. In this case we easily get 
the relation 

[L,L] = 1. (6.14.7) 

The string equations in the form of Takasaki allows us to notice the similarity to the different 
problem. The dispersionless limit of ( |6.14.7| ) (and also of ( |6.14.5| ), (|6.14.4| ), where r(M) = M^, 
and of (|6.14.6| )) will be written as 



AA 



n e Z, 



(6.14.8) 

/i = /i. (6.14.9) 

The case when A and p, are complex conjugate of A and of fj, respectively, is of interest. These 
string equations (mainly the case n = 1) were recently investigated to solve the so-called 
Laplacian growth problem, see |H3|. We are grateful to A.Zabrodin for the discussion on this 



problem. For the dispersionless limit of the KP and TL hierarchies see [31 . 
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In case the function r{n) has zeroes (described by divisor m = (Mi, ...): r{Mk) = 0), one 
needs to produce the replacement: 

A^A(m), A^A(m), (6.14.10) 

where new matrices A(m),A(m) are defined as 

(A(m))„- = 5,,,„i, J ^ Mfc, (A(m)),,, = 0, j = M,, (6.14.11) 

(A(m)),j = <5,,,+i, z ^ Mfc, (A(m)),j =0,i = M^. (6.14.12) 
This modification describes the open TL equation (??): 

dtA*,Vn = 5(n)e'^"-i"'^" - 6{n + l)e'^"-'^"+i. (6.14.13) 

The set of fields (^oo, • • • , ^a/i, 0a/i+i • • • consists of the following parts due to the conditions 

Ms Mk+i A/i+1 

5] <^„ = 0, ^ 0„ = 0, E 0„ = 0, (6.14.14) 

n=-oo n=A4'fe+l n=oo 

which result from r(Mfc, t,t*) = 1. 

Remark 2 T/ie matrix r(M) contains (M — 6j) m i/ie denominator. The matrix (M — fei)^"*^ 
is i^(A — bi)^^K^^{l + 0(t)) (compare the consideration of the inverse operators with j^/j. 



The KP tau-function (|3.4.8|) can be obtained as follows. 

G(M,t,T,t*) = G(0,T,0)f/(M,t,r), f/(M,t,t*) = f/+(t)f/-(M, t*). (6.14.15) 

f/+(t) = exp(e(t,A)), f/-(M,t*) = exp(e(t*,Ar(A + M))) , (6.14.16) 

The matrix G{0,T,0) is related to the transformation of the eg. ( |6.3.1^ ) to the eg. ( |6.3.lTD . 

By taking the projection U ^ U for non positive values of matrix indices we obtain a 

determinant representation of the tau-function ( ^.4.8} ): 

detf/__(M,t,t*) . *x , . 

r, M,t,t* = s r = detf/„_ M,t,t* , 6.14.17 

det([/^_(t))det(f/r_(M,t*)' V ' ' ^ I ) 



since both determinants in the denominator are egual to one. Formula (|6.14.17|) is also a Segal- 
Wilson formula for GL(oo) 2-cocycle Cm iU^{—t),U~{—t*)). Choosing the function r as 
in Section 3.2 we obtain hypergeometric functions listed in the Introduction. 



Remark 3 Therefore the hypergeometric functions which were considered above have the mean- 
ing of GL{oo) two-cocycle on the two multi-parametric group elements U~^{t) and U~ {M,t*) . 
Both elements U^(t) and U^{M, t*) can be considered as elements of group of pseudo-differential 
operators on the circle. The corresponding Lie algebras consist of the multiplication operators 
{z^; n G Nq} and of the pseudo-differential operators {(^jf{z-j^ + M)^ ; n G A'^o}. Two sets of 
group times t and t* play the role of indeterminates of the hypergeometric functions ( 6.12.^ ). 
Formulas ( ^.4-8{ ) and ( \6. 3. 8^ ) mean the expansion of GL{oo) group 2-cocycle in terms of corre- 
sponding Lie algebra 2-cocycle 

CM{z,-r{D)) = r{M), CM{r{D)z,-r{D + M)) = r{M)r{M). (6.14.18) 

2/ Z 
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Japanese cocycle is cohomological to Khesin-Kravchenko cocycle /j37|/ for the ^DO on the circle: 

Cm ~ Co ~ cum, (6.14.19) 

which is 

uom{A, B) = ^= / resQA[\og{D + M), B]dz, A,B e "^DO. (6.14.20) 

2'Ky/ — l J 

For the group cocycle we have 

CM(e-i:^"*",e^i:(^"'^'(^))"*") =r,(M,t,r), (6.14.21) 
where we imply that the order of^DO r{D) is 1 or less. 



Remark 4 It is interesting to note that in case of hypergeometric functions pFg ( |<§. 0.1'^ the 
order of r is p — s (see Example ?>), and the condition p — s < 1 is the condition of the 
convergence of this hypergeometric series, see / |^/ . Namely the radius of convergence is finite 
in case p — s = 1, it is infinite when p — s < 1 and it is zero forp — s > 1 (this is true for the 
case when no one of in ( \8. 0. 1 1\ ) is nonnegative integer). 



Remark 5 The set of functions {w{n., z),n = M, M + 1, M + 2, . . .}, where 

w{n,z) = exp (- f; C (-r{D)Y] ■ z\ (6.14.22) 



m=0 "-^ 



may be identified with Sato Grassmannian l \6.13.d{ ) related to the cocycle l \6.14.21\ ). The dual 



Grassmannian \6.13.(\ ) is the set of one forms {ty*(n, 2;), n = M, M + 1, M + 2, . . 



w [n,z) = explf^r^f-ri-D)) ] ■ z~'^dz. (6.14.23) 



\m=0 



Z 



7 Appendices B 

7.1 A scalar product in the theory of symmetric functions 

We consider a ring of symmetric functions of variables Xi,i = 1,2, ... ,N , where N may be 
infinity. Let us remember the notion of scalar product. Different choices of scalar product give 
rise to different systems of bi-orthogonal symmetric polynomial functions (like Schur functions, 
Macdonald functions and so on). 

Actually we consider two sets of variables, the first one is x = (xi, 0:2, . . . , xat), the second 
is p = {pi,P2, ■ ■ ■)■ These two sets are related via 

N 

Pm = T.^T, rn = 1,2,... (7.1.1) 

i=l 

If one takes 

eE:.i ^-"IW = y: Px{p)Qx{p*) = E (7-1-2) 
then for any partitions /i, A 

{P^.,Qx)v = 5^,x (7.1.3) 
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with respect to the scalar product 

{p,,Px)v = z'^x^5,,x (7.1.4) 

When V = 1 then 

z?=zx, zx = l[^"''m,\ (7.1.5) 

i 

where rrii = mi{X) is the number of parts of A equal to i 

Pt^ = I[Pl^^^ Pt^, = J2^k\ = Pm (7.1.6) 

i k 

If one takes 

= Yzr^ (7.1.7) 

he gets Macdonald's polynomials Qa(p) = Qa(p; Q', i)^ Pxiv) = Px{P'^ li t)- this case 

z'^^^ =z,{q,t) = z,Y[—-^, z, = Y[t"^'m^ (7.1.8) 

i=l i 

where /(A) is the length of A and = mj(A) is the number of parts of A equal to i. 

7.2 A deformation of the scalar product and series of hypergeomet- 
ric type 

For different choice of function v we choose numbers a^'"\b^'"^ and a function r^""^ which is a 
function of one variable M, say. We define 



(M) = n r^^KM + a^'^^z + b^^^j) (7.2.1) 



as a product over all boxes of the Young diagram. In what follows we shall omit the superscript 
(v). 

Now let us deform the scalar product (|7.1.3|) in a following way 

{P>.,Qx)v,r,M = rxiM)6^,x (7.2.2) 
To get series of hj^ergeometric type we consider 

(e^™=i^fc=i fc m ^ e^'"=i " ™ /D.r.M (7.2.3) 

Using (|7T1) 

(eE:=iEr=i''- -r^,eS™=iS^=i''™ ^^"''^),,,,M = J2rx{M)Px{^'')Qx{y'') (7.2.4) 

A 

In this paper we shall consider two cases of ( [7. 1 . 7| ) . First one is 

Vn = 1 (7.2.5) 

for all n. The second one is 

V2n = 0, V2n+1 = 2 (7.2.6) 

The first case corresponds to Schur function series. The second one corresponds to projective 
Schur function series. These are the cases when the scalar product (|7.2.2| ) can be presented as 



a vacuum expectation of certain fermionic fields, see Section 4 and ||52| 
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7.3 Scalar product of tau functions. A conjecture. 

Before we considered tau functions of hypergeometric type (|3.4.8| ). In this subsection we shall 
consider more general TL tau functions. Let ri(t, t*), r2(t, t*) be TL tau functions of the form 

ri(t,t*) = e^(*)-eS™=i'"*'"*™, r2(t, t* ) = e^^**) ■ eS-=i (7.3.1) 

where e^ = l + X + -- -,e^ = l + Y-\----, and where the operators X{t) and Y{t*) are linear 
combinations of vertex operators (|6.1.1|) and ( |6.1.2|) respectively: 

X,(t) = J f{z,z')VUt,z)V^{t,z')dzdz' (7.3.2) 
y,(r ) = J g{z, z')Vo{t*, z)V^{t*, z')dzdz' (7.3.3) 



If the TL tau functions (^XTJ) have the typical forms ([LaOD (see , J29[ , [0 ) 

ri(t,r) = Y.Kx^sx{t)s^{t*), r2(t,t*) = Y.My^^sx{t)s^{V) (7.3.4) 

A, /J A,/i 

then their (standard) scalar product ( |2. 1.151 ) is obviously equal to 

< ri(t,7),r2(7,t*) >= ^ iVA^SA(t)s^(t*), iV^^ = J^Ka.M,^ (7.3.5) 

It is not immediately clear that what we got is a tau function again (without solving Hirota 
equations). However on the other hand we have 

< ri(t,7),r2(7,r) >= e^^*) ■ e^^**)- < eS-=i eS™=i >= e^^*) ■ e^^**) ■ eS-=i 

(7.3.6) 

The last relation shows that we get tau function again according to general results of |^. (It 
was shown there that the action of vertex operators e-^ to KP tau function gives rise to a new 
KP tau function. TL tau function is the tau function of the pair of KP hierarchies related to 
the sets of times t and t*, see |l27[|). 

Due to ( p.7.1| ) it means that if the following integral over complex planes 7^, = 1, 2, . . . 



/ n(t,7)e-S™^.'"l-lV2(-7,t*) n ^ (7-3.7) 

•' m=l ^ 

exists, then it is a new TL tau function. Here the variables 7 = (71,71,...) are complex 
conjugated to 7 = (71, 72, . . .). 

The proof can be done also directly as follows 

/ ri(t,7)e-i::^.H7.lV^(_^,t*) n ^ = (7.3.8) 

m=l ^ 

gX(t) . gy(t*) . y gX;"=i'"*-Tm-m|7„|2-m7„f^ -Q md^lm ^ (7.3.9) 

m=l ^ 

^(t) .e^(**) .6^^=!'"*'"*™ (7.3.10) 



m=l 

which is tau function. 
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Let us try to consider more general situation. Let ro(t,t*) be a TL tau function with 
the following properties. First, ro(t,t*) = ro(t*,t). Second, each integral over complex plane 
7m,m = 1,2, . . . 

j Toil, -l)d^lni = dm (7.3.11) 

where 7^ is complex conjugated of 7^, is convergent and non- vanishing. Then 

m=l 

Given tq, we introduce the following scalar product 

m=l 



/-o(7,-7)n^ = l (T.3.12) 



< f,9 >o= / /(7)ro(7,-7k(-7) IT — ^ (7-3.13) 



Let us present a conjecture that scalar product of TL tau functions is TL tau function again. 
Conjecture. Let ri(t, t*), T2{t, t*) be TL tau functions. If the integral over complex planes 
7m,m = 1,2, . . . 

/ ri(t,7)ro(7,-7)r2(-7,t*) U (7.3.14) 

m=l ^rn 

exists, then it is a new TL tau function. □ 

7.4 Integral representation of scalar products 

For the scalar product ( |7.2.2| ) one can present the following integral realization. 
In case r = 1 

m=l '^iTT- 

where mtm are power sum functions, mim are complex conjugated variables, and integration is 
going over complex planes {tm,im) ■ Normalization factors am depend on m: 

am = - (7.4.2) 
m 

In case r(n) = n there is a representation 

U,9)r.n =[■■■[ fiz)g{z)e-^l^\^^\'\Aiz)\' J] ^ (7.4.3) 

where z = {zi, . . . , z^) are complex conjugated to z = {zi, . . . , z^), normalizing constants are 

bk = rr (7.4.4) 

A is Vandermond determinant 

N 

A{z) = liizk-Zi) (7.4.5) 

i<k 

and one puts A.{z) = 1 for = 1. The integration is going over complex planes {z^, z^). 
For the scalar product ( |7.1.4| ) we present the representation 



{f,g%= [■■■ [ /(%(t)e-S:-e*-*- n (7-4.6) 
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where mtm are power sum functions, mtm are complex conjugated variables, and integration is 
going over complex planes (tm,^m) • Normalization factors depend on m: 



2ttv„ 



m 



In particular for Macdonald's polynomials Qx{p; q, t), Pa(p; Q, t), see (|7.1.7| ), we get 



1-9 



Y^oo _ 1- 

Qx{p; q, t)PxiP; g, t)e~ ^--i p™p™;^ 

1 - 



(T^ H dPradpm ^ ^^^^ 



m=l ""m, 



am. = 271- 



mil-t" 



(7.4.7) 

(7.4.8) 
(7.4.9) 



8 Appendix C. Hypergeometric functions 

Ordinary hypergeometric functions First let us remember that generalized hypergeometric 
function of one variable x is defined as 



pF, (ai, . . . , Op; 6i, . . . , 6,; x) = }^ — TT^-y 

n=0 \"i)n ' ' ' K'JsJn 



(8.0.10) 



Here (a)„ is Pochhammer's symbol: 

T(a + n) 



(«)r 



r a 



a(a + 11 



(a + n — 1). 



(8.0.11) 



This hypergeometric function solves the equation 



^ n ( + '^i ) I p-^s ('^1' • • • , Op; • • • , 2;) 



0, 



(8.0.12) 

Given number q, \q\ < 1, the so-called basic hypergeometric series of one variable is defined 



as 



, , , ^{q''';q)n---{q''^;q)n 

bi,...,bs]q,x) = 2.7-^ 



^0 i.(f^'^(l)n---{q^']q)n {q;q)n' 
Here (g°, g)„ is g-deformed Pochhammer's symbol: 

(6; g)o = 1, (6; g)„ = (1 - 6)(1 - 6g^) ■ ■ ■ (1 - bq^-'). 



(8.0.13) 



(8.0.14) 



Both series converge for all x in case p < s + 1. In case p = s + 1 they converge for \x\ < 1. 
We refer these well-known hypergeometric functions as ordinary hypergeometric functions. 



The multiple hypergeometric series related to Schur polynomials ||26|| , ||25|] , [p8 



There are several well-known different multi- variable generalizations of hypergeometric series 



of one variable If one replaces the sum over to a sum over partitions A = (rii, n2, . . .), 

and replaces the single variable x to a Hermitian matrix X, he get one of generalizations of 
hypergeometric series which is called hypergeometric function of matrix argument X with indices 
a andh m],^: 



/ cti, . . . , Op 



X 



E 



{ai)x ■ ■ ■ iap)x Zx(X.) 



. (&i)a---(Ma |A|! 

i(A)<JV 



(8.0.15) 
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Here the sum is over all different partitions A = {rii, n2, . . . , n^), where ni > n2 > ■ ■ ■ > Uk, 
k < |A|, |A| = rii + ■ ■ ■ + Ur and whose length /(A) = k < N (n^ 7^ 0). X is a Hermitian N x N 
matrix, and Zx(X.) are zonal spherical polynomial for the symmetric spaces of the following 
types: GL{N, R)/ SO{N), GL{N,C)/U{N), or GL{N, H)/ Sp{N) see [||. The definition of 
symbol {a)\ depends on the choice of the symmetric space: 

1 k — 1) 
(a)A = (a)ni(a )„2 ■ ■ ■ (a )nk, (a)o = 1, (8.0.16) 

a a 

where a = 2, a = 1 and a = | for the symmetric spaces GL{N, R)/ SO{N), GL{N,C)/U{N), 
and GL{N, H)/Sp{N) respectively. The function ( |8.0.15| ) actually depends on the eigenvalues 
of matrix X which are = [xi, . . . ,xn) . In what follows we consider only the case of 
GL{N,C)/U{N) symmetric space. Then zonal spherical polynomial Zx(X.) is proportional to 
the Schur function sx{xi,X2, --.yXN) corresponding to a partition A Q, s\{xi,X2, ■■■,xn) is a 
symmetric function of variables Xk- 

For this choice of the symmetric space the hypergeometric function can be rewritten as 
follows 



s 



Qiii ■ ■ ■ 1 dp 

bi,...,bs 



A 

i(A)<iV 



where Hx is 'hook product': 

i^A = n ^ii' = i^i + n'j-i-j + l) (8.0.18) 

and 

(a)A = (a)„^(a - 1)„2 ■ ■ ■ (a - + 1)„^, (a)o = 1 (8.0.19) 

Taking = 1 we get (|8Xl0|) . 

For this hypergeometric functions we suggested different representations like ( |6.11.12| ),( |6.11.13D 
or like (|6.12.14|) , (|6.12.15|) , (|6A21^ ). 



Let |g| < 1. The multiple basic hypergeometric series related to Schur polynomials were 
suggested by L.G.Macdonald and studied by S.Milne These series are as follows 



(ai,...,ap;6i,...,6,;g,x = }^ — — ^ — — — — — sa (x , 8.0.20 



\ 

l(X)<N 



where the sum is over all different partitions A = {rii, n2, . . . , n^), where ni > n2 > ■ ■ ■ > rik > 0, 
k < |A|, |A| = rii + ■ ■ ■ + Hk and whose length /(A) = k < N. Schur polynomial sx (^^)) with 
N > /(A), is a symmetric function of variables x^ and defined as follows 0]: 

s,(x^) = ^A = det(a;rOi<„<7v, 5 = (iV - 1, iV - 2, . . . , 1, 0). (8.0.21) 

as 

Coefficient (g'^; q)x associated with partition A is expressed in terms of the g-deformed Pochham- 
mer's Symbols {q^;q)n (|8.0.14| ): 

{q"; q)x = iq'; q)nAq'-\ qU " " " q)n,. (8.0.22) 

The multiple g"*^'^) defined on the partition A: 

= (8.0.23) 
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and g-deformed 'hook polynomial' H\{q) is 



{*,i)eA 



hi 



(8.0.24) 



where A' is the conjugated partition (for the definition see 0). 

For this hypergeo metric functions we suggested different representations like ( |6.11.14| ),( |6.11.l5D 
or like (|6.12.29|) , (|6.12.3q) . 

For = 1 we get (|XT§. 

Let us note that in the limit g — > 1 series ( |8.(J.2U| ) reduces to ( |8.U.17| ), see 
Hypergeometric series of double set of arguments ||26|| , []25|] , ^8 



Another generalization of hypergeometric series is so-called hypergeometric function of two 
matrix arguments X, Y with indices a and b: 



oJF g (oi, . . . , Op; 6i, . . . , bg] X, Y) — ^ 



{ai)x---{ap)xZx{X)Zx(Y) 
(&i)a---(&s)a |A|!Z,(I^) • 



(8.0.25) 



Here X,Y are Hermitian N x N matrices and Zx(X.), Zx{Y) are zonal spherical polynomials 
for the symmetric spaces GL{N,C)/U{N), GL{N, R)/SO{N) and GL{N,H)/Sp{N) see fH. 
The notations are the same as in previous subsection. Again we shall consider only the case of 
GL{N, G)/U{N) symmetric spaces. In this case (|8.0.25| ) it may be written as 

{ai)x ■ ■ ■ iap)x SA(x^)sA(y^) 



pj^s (ai, . . . , ap; 6i, . . . , bs] x^^^ y 



[N) 



E 



(6i)a---(6s)a Hxsxiin 



(8.0.26) 



where 1^ = (1, . . . , 1). 

We suggest a different representations like ( |6.12.2U| ) or different integral representations. 
The q-deformation of the hypergeometric function (p.0.26| ) is as follows 











■ ,bs 


g,x^,y^) 



(g"^ g)A • • • {q""; q)x g"^^^ (x^) ^a (y^) 

. {q'^; g)A ■ ■ ■ {q'-,q)x Hx{q) Sx (1, g, q^ g^'^) 

;{A)<jv 



E 



(8.0.27) 



This is the multi- variable basic hypergeometric function of two sets of variables which was also 
studied by S.Milne, see p5[ . 



There are also hypergeometric functions related to Jack polynomials C^x"' IP5|| : 





[d) 


(ai, . 


. .,ap;6i,...6^;x^,y^) = 


E 

A 


(ai 


^{d) 
l\ ■ ■ 


.(ap)f Cf(x^)Cf(y^) 




^{d) 
)\ ■ ■ 


■(6,)f |A|!Cf (P) 



where 



Here (c) 



liX) 

n 

i=l 



(8.0.28) 



(8.0.29) 



rii 



c(c+l) 



{c + k — 1). For the special value d = 2 the last expression (|8.0.28|) 
coincides with (|8.0.20| ), and reduces to (|8.0.26|) as |g| — > 1. The open problem is to get a 
fermionic representation of (|8.0.28|) for arbitrary value of the parameter d. 
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9 Further generalization. Examples of Gelfand-Graev 
hypergeometric functions 

9.1 Scalar product and series in plane partitions 

Let us consider the scalar product (|3.3.9|) 

{Sf„sx)r,n = rx{n)Sf,^x (9.1.1) 

For the simphcity of notations sometimes we shall write (, )r instead of (, )r,n- 
We introduce notations 

oo oo oo 

(t>P) = J2 P^^^' (P' A) = J2 Pmf^ini, (P, A*) = J2 Pmf^im (9.1.2) 
m=l m=l m=l 



Lemma 



(SA, e(P'^)s^),,„ = (s^e^P'^), sx) = rx{n)sx/M (9.1.3) 



One sees that for /? = and r = 1 the right hand side is equal to 5^^x- ^ 
Proof One develops e'-^'^*-' using 



e 



E^-(7)^.(A) (9.1.4) 



and taking into account the orthogonality of Schur functions ( |9.1.1|) and the formulas (see 
section 5 of chapter I of Q]) 

SvS^, = (^tuSx, sa/^ = ^^u^^ (9.1.5) 
A y 

□ 

Let us consider a set of scalar products ( |9.1.1| ) which differ by the choice of r. The scalar 
product of functions of Pi will be labelled by script i: we shall denote this scalar product as 

Lemma 

e{P.,ft+7.-i) = ^ ,,^(7,)s,,/,,(A)s.,(7.-i) (9.1.6) 
and for given partition z/j_i we have 

where 

7. = (^,f,f,...) (9.1.8) 
□ Proof follows from the relations (see Q]) 

(9.1.9) 

and definition ( p.3.9| ). □ 

Let us consider a set of {pj, z = 0, 1, 2, . . . , /c + 1} and a set of i = 1, 2, . . . , + 1} with 
conditions 

Po = 0, Pfc+i = p (9.1.10) 
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With the help of Lemma we consider scalar product of scalar products (compare with 

mi))-- 



(e(P.A+7.-i) . . . ^ ^g(P3,/33+72)^ ^^(p„;3,+,0^ e(Pi,/3i+7o))^^^^^^^)^^^^^^^ . . .)^,,_,,^^^ (9.1.11) 

= E ^^t:!n^.~i)r^t?K-2)---^l?K)r«(nO 

(/52)s.,/.o(A)s.o(7o) (9.1.12) 
= 4(p,r^-\A'^-\/3,Po) (9.1.13) 

where r^^^^ = (r(i),r(2), . . . , r^''^^)), Afc_i = (rii, . . .,nk-i), 5^,0(70 = 0) = 5o,i^o- 
This is the sum over all plane partitions ^ with the largest part k: vr^'^) = (z^i ^ ■ ■ ■ C u^). 
Now let us consider scalar product of different / and with respect to standard scalar product 



(|9.1.1| ) . Putting po = we get a function oi (3 = {/?m, 2 = 1, . . . , A;, n = 1, 2, ...},/?* = 



1, . . . , /, n = 1, 2, . . .}, which will appear to be a tau function 

rf,r(A;_i,A*„/3,/5*) = (4(p,ffe_i,A0),/,(p,r;_i,/3*,0)),w,„^ (9.1.14) 

where f = (f(-'-\ f(^), f^^" -*-)), r = (r^^), r^^^ . . . , r^'^"-'-^ r'^'^^) and A;_i = (ni, . . . , n;_i), A^ = 
(n^, . . . , n^). We see that ( |9.1.14| ) is a sum over a pair of plane partitions. 



9.2 Integral representation for r 

In case the function r has zero the scalar product (|9.1.1|) is degenerate one. For simplicity let 
us take r(0) = 0, r{k) ^ 0, k < n. The scalar product ( |9.1.1D is non generate on the subspace of 
symmetric functions spanned by Schur functions {sa,/(A) < n}, Z(A) is the length of partition 
A. 

(/, 9)r,n = / • ■ ■ / n dzidzifir{ziZi)Aniz)Aniz)f{z^)g{z^) (9.2.1) 

i=l 



where z"' = {zi, . . . , Zn), z"- = {zi, . . . , Zn) and 



An{z) = l[{zi - %), A„.{z) = llizi - zj) (9.2.2) 

i<j i<j 



Thus get the following integral representation for ( |9.1.14|) 



rf,r(Ai-i,A^,A/5*)= (9.2.3) 

l-l "j k " 

Yl n dzjidzjififU){zjiZji)DjS{zki - Zi_-i^^i)S{zki - ^f_i_i) 

j=l i=l j=l i=l 



j ■■■ j n n dzjidzji^fU){zjiZji)DjS{zM - z^_-^^i)S{zk^ - f[ Y[ dz*idz*ifi^u){z*iZ*i)D* 

j=l i=l j=l i=l 

(9.2.4) 

D, = A„^.(4)A„^.(4)e«^ D* = A„j(z;)A„^(4)e«^* (9.2.5) 

oo / 1 "-j-l \ CO ni 

= E EM'" + - EI^.-m)"" '^' > 1' ^1 = E E(^ifcrA- (9-2-6) 



m=l k=l \ ^ i=l / m=l k=l 

1 



= E E /^In. + - E i^UT , J > 1, = E E(^tfc)'"/5L (9.2.7) 



m=lfc=l \ "^4=1 / m=lfc=l 
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9.3 Multi-matrix integrals 



In case A^\n) = n + aj,r^^\n) = n + a* we get that the tau function ( 9.2.3 ) is equal to the 
following multi matrix integral 



■ ■ I n dM,dMje^''^^'^^e^^6{Mk - M;_^)6{Mk^ - M;_i,J f[ dM*dM*e'''^^^h^^ (9.3.1) 

OO / -1 N oo 

V, = Tr{M,r (Pjm + -Tr{M,.^r) , j > 1, 1^1=^: Tr(Mi)"^A^ (9.3.2) 
ni=i \ m J 

oo / -1 x oo 

V;=Y. TriM;r [P^m + -Tr{M*_,r) , j > 1, V* = ^ Tr(M;)"^/5L (9.3.3) 



TTl 

m=l m=l 

where Mj, Mj has a size rij = M + aj and M*, M* has a size n* = M + a*. 



9.4 Fermionic representation [[T2| 



Due to (13.3.101 ) scalar product (|9.1.14|) is equal to a vacuum expectation value. Let us write 



down it explicitly using |]12[ . 

Formula (|6.3.1|) is related to 'Gauss decomposition' of operators inside vacuums (M| . . . \M) 
into diagonal operator e^"'-^^ and upper triangular operator e^*^*^ and lower triangular operator 
g-j?*(f ) ^Yi^Q last two have the Toeplitz form. Now let us consider more general two-dimensional 
To da chain tau function 

r = (M|e^(*)^e-^(**)|M), (9.4.1) 
where we decompose g in the following way: 

(7(/?,/5*) =e^i(^i)-- •e^^(''^)e-^'(^')-- -6-^1(^1*), (9.4.2) 

where 

/3* = (Ai,A2,...), (3: = i/3:,,(3:,,...) (9.4.3) 

00 00 

i,(A) = AtmPkm, AkiP\) = E Akn.P\m (9-4.4) 

m=l m=l 

Here each of Ak{P*) has a form as in ( ^.5.3] ) and corresponds to operator r^'^^D), while each 
of Ak{/3*i^) has a form of (|3.5.4|) and corresponds to operator f^^\D): 



A 



Akm ^ 



\ r m 
— ■= ii)*{z)(f^''\D)zY i){z), m = l,2, ... (9.4.5) 

1 /■ /I \ ™ 

7/^*(2) (^-r('=)(D)j V'W, m = l,2,..., (9.4.6) 

Collections of variables j3 = = play the role of coordinates for some wide 

enough class of Clifford group elements g. This tau function is related to rather involved 
generalization of the hypergeometric functions we considered above. Tau function ( p.4.1| ),( p.4.2| ) 
may be considered as the result of applying of the additional symmetries to the vacuum tau 
function, which is 1, see Appendix "The vertex operator action". 

Let us calculate this tau function. First of all we introduce a set consisting of m + 1 partitions: 
(A\...,A™,A'"+^ = A), = A° C A^ C A^ C ... C A™ C A'"+^ = A, (9.4.7) 
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see [Q for the notation C for the partitions. The corresponding set 

= (A\0\...,^™), e' = y^^-x\ i = i,...,m (9.4.8) 

depends on the partition A and the number m + 1 of the partitions. We take as se(t*,P*) 
the product which is relevant to the set 9^ and depending on the set of variables fii = {fJ^ij} 
(^ = (l,...,m + l),j = (l,2...)) 

se'^if^) = ■SAi(/Ui)sei(/U2) ■ ■ ■ serr^{fim+i)- (9.4.9) 
Here s^i is a skew Schur function (see Q). Further we define function r0m(M): 

re^{M) = r^i{M)rl^{M) ■ --r^AM), (9.4.10) 
where the function rgi(M) , a skew analogy of rx{M) of (|2.2.1|) , is 



s 

rg,{M) = l[r{nf - j + 1 + M) ■ ■ ■ r{n^^^^^ - j + M), (9.4.11) 
i=i 

where A*"*"^ = {ni~^^\ . . . , n^*"*"^)). If the function r^{m) has no poles and zeroes at integer points 
then the relation 

-MM) = ^|tii^, ^ = l,...,m (9.4.12) 
is correct. To calculate the tau function we need the Lemma 

Lemma 3 Let partitions A = {ii, . . . ,is\ji — 1, . . . ,js — 1) and A = {ii, . . . ,ir\ji — 1, . . . ,jr — 1) 
satisfy the relation A ^ A. The following is valid: 

(OIV^J ■ ■ ■ V'-.v ■ ■ ■ ^-J,e^"(^'V-,, ■ ■ ■ ■ ■ ■ |0) = 

= (_l)ii+-+i^+ii+-+i.sg(;5.)^g(0), ^ = A-A. (9.4.13) 

(OIV^* ■ ■ ■ i^li^-,. ■ ■ ■ ^-,,e^'(^?)^l.^ ■ ■ ■ ■ ■ ■ ^ijO) = 

= (-l)j'i+-+>+Ji+-+J^^se(/5i)re(0), ^ = A - A. (9.4.14) 

□ Proof: One proof is achieved by a development of,e"^ = l + A + -- -,e"^ = l + yl + -- - and a 
direct evaluation of vacuum expectations ( |9.4.13D ,( |9.4.14| ). (see Example 22 in Sec 5 of |^ for 
help). 

The second proof of (|9.4.13| ),( |9.4.14|) is achieved using the developments (|3.5.17| ), (|3.5.18 



respectively. □ 

Then we obtain the generalization of Proposition 1: 

Proposition 18 



rM(t, t*; /?,/?*) = E E E feiWr^i (M).ej(t, P)s^.^{t\n, (9-4.15) 
where Tqu^M) and Tqi (M) are given by ^9.4-l^ )- 
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With the help of this series one can obtain different hypergeometric functions. 
In the end of the subsection we put j3 = 0. For the case t = t(x^), all /3 = one can obtain 
the analog of 

(9.4.16) 



where 



+ 00 



D = z-^, r'k{D) = n{-D) 



(9.4.17) 
(9.4.18) 



In ( p.4.16| ),( |9.4.17] ) are operators which act on z variable of 

In cases t = t(x^) (using (|X|) and (|3A^ )) and t = -t(x^) (using (|6Xl0| ) and (|3X^ )) 
one gets the following representations 

Proposition 19 

nj<7 i-^i -^i) 



t{M, t(x^), t*; (3 = 0, 13*) = A-^e^' ■ ■ ■ e''^ A, A = 
r(M, -t(x^), t*; (3 = 0, (3*) = A'^-"^' ■ ■ ■ e-^'e-'^A, A 



(xi ■■■xn) 



N-l+M 



where 



N N 

V = Y.^r{t\Xi), r/fc = 5Z^rfe(t*,Xi) 
i=l i=l 
N N 

V' = J2^r'{t*,Xi), r]k = ^C,.,(t*,Xi) 



i=l 



i=l 



and (of ^TD) 



+ 00 



+ 00 



_d_ 

dxi 



(9.4.19) 
(9.4.20) 

(9.4.21) 
(9.4.22) 

(9.4.23) 



{r',(t*,Xj) = ^ (xirfc(A))'", irSt*,Xi) = tm{xirk{-D,0)"' , Di = Xi 

m=l m=l 

Let US note that [^^'^(t*, Xj), ^r'fe(t*, x^)] = and [^rfe(t*, a;,), ^rfc(t*, Xj)] = for all k,i,j, while 

in g eneral [r]k,r]n] 7^ 0- 

( 19.4.201) may be also obtained from (|9.4.19|) with the help of ( |3.4.10|) . 

Looking at (|9.4.19|) one easily derives a set of linear equations, which are differential equa- 
tions with respect to variables (3* and t*: 



_ e'?' ... e'"= 



■ [J2 ^^rkm\ ■ e-'"=-' ■ ■ ■ e-^'' j (Ar(M, t(x^), t*; P = 0,P* 
1^ - e''' ■ ■ ■ e'^' ■ (E ^^rk{D,)^ ■ e-^^' ■ ■ ■ e-^'' j (Ar(M, t(x^), t*; P = 0, (3*] 



^ ^ ^ 

km \i=l / / 

(9.4.24) 

= 

(9.4.25) 

where A is the same as in ( |9.4.19| ). 

Proposition 20 For the tau function Tr{M, t, t*) of l \3.4-(^ and the vertex operators 
defined by ( \6. d.ij ) we have 

e-^r(ft) . . .e-^,-'(^^) . . ■e^'^'^'^'*) ■ ■ -e^'^'^^^*) ■ r,(M, t,t*) = r(M, t, t*; /3, /3*), (9.4.26) 

. . .e^:°'(/3r) . ..^-^f^(M . . ^^(M,t,t*) = r(M, t, t*; /?, /?*). (9.4.27) 
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9.5 Loop matrix integrals 

Let us consider three matrix integral used in BTI 



J = J gE:.i(/^"^Trx::>+/5;;,Trx™)^V^Trx_x+g-,^Trx_nx+f^-^^_^^^^^ ^g_5_^) 



Here X± are Hermitian n by n matrices, and Q G U{n). (In ( |9.5.1| ) we use jSm (and (3^) instead 
of tm (and of in pTj). 

Using ( pXTTD , (pX8| ) and ( pXSD we get 



r_ V- „ ^^^ < gA(X_)g^(X_),g^(X+)g^(X+) >,,„ ^ _ ^ {n)xSx/u{f3)sx/u{P*) , . 



As we see / =<< r >> (see ( p.6.8|) ), where 



where 

^ = (_g)«e^('^)e^(^*)e- '"'^(""^ S™=— '"^'^'"'^'"^ (9.5.4) 

where A is the same as in the tau function for the two-matrix model (r(m) = m). We have 
I = TrF„^ = Ea^aa- 

9.6 The example of Gelfand,Graev and Retakh hypergeometric se- 

m 



ries 



Below we also put j3 = 0. Let us consider the tau function: 

r(M, /3, /?; 0, = {Mle^^^^e-"^'^^'^ ■ ■ ■ e-^'^^'^^e-^^^^\M). (9.6.1) 

We put 

2 3 

/3=(a:,y,y,...), /? = (l/i, 0, 0, . . .), P* = {y,+,,0,0, . . .) t = {l,...,l). (9.6.2) 
We obtain the series 

T{M,x,y,,...,m^,) = 2 r(„,^...+„,^,)(M)re. (M) " " " (^^'+;)"'"^ = (9.6.3) 

ni,..Xi=o rii!---n,+i! 

J2 c m, . . . , riz+i xyi "1 ■ ■ ■ Or2/i+ir'+\ c{ni, . . . , = — — — ^—,(9.6.4) 

ni,...iCi6^ r(ni + l)---r(r2i+i + l) 

where 0';)^ corresponds to the set of simple partitions-rows 

= (ni), = (ni + n2), . . . , A/+i = (ni H h n^+i) (9.6.5) 

When functions 6j(ni, . . . , m+i) defined as 

6j(rai, . . . ,n/+i) = r , z = !,...,/ + ! (9.6.6) 

c{ni, . . .,ni+i) 



68 



are rational functions of (ni, . . . , n/+i), then tau function ( |9.6.3D is a Horn hypergeometric series 



25 



Above series for the special choice of functions r' (D) can be deduced from the Gelfand, Graev 
and Retakh series |3 defined on the special lattice and corresponding to the special set of 
parameters. Let us take the rational functions r^{D): 



n(») / ^ . (^)^ 



rUD) = ^^=^^^^"{./ , (t = 0,...,l), r\D)=r(D) (9.6.7) 



Let define = + s(°) + 2 E5=i(p^^^ + 2s(^')) + + 5^'+^) + / + 1 and consider complex 
space C^. In this space we consider the / + 1-dimensional basis B and the vector v consisting 
of parameters. 

po = so = 0, p,=p(^-i) + s», s, = s(^-^)+p«, ^ = (1,...,/) 

i+i 

Pi+i = +/+^\ = s« + AT = Y^ip^ + sj) + l + l (9.6.9) 

i=i 

f ~ ~(^Po+soH hPi-i+Si_i+l + ■ ■ ■ + Gp^+siH |-Pi_i+Si_i+Pi) + 

+ {^Pi+si+-+pi.i+Si.i+pi+i + ■ ■ ■ + ep-^j^.sj^+-+pi_i+si_i+pi+si), ^ = 1, ...,/ + 1 (9.6.10) 

i 

where ej = (0, . . . , 0, 1, 0, . . .). The lattice B G is generated by the vector basis of dimension 
^ ^ ^ 

N 

y = f + + +eN-i-i+u « = !,...,/ + ! (9.6.11) 



Vector V G is defined as follows (compare with ( |9.6.10|) ): 



i / («— 1) I I (*— 1) I 

~ ^PQ+soA hPi-i+Si-1+1 + ■ ■ ■ + a'p{i-i)Gpg+soA hPi-i+Si-i+p'*-^) 

+ '^1 6po+so + --+Pi-l+Si-i+p('-^'+l + ■ ■ ■ + C'<j(i)^PO+so + ---+Pi-l+Si-l+pJ + 

+ ((^1 ^ - l)ep(,+so+-+p,-i+s,_i+p,+l H ^ ^ l)6po+so+-+p>-l+s>-l+p,+s('-l) + 

+ — l)ep|j_,_g|j_,_...^_p._^_,_^._^_,_p^_,_g(i-i)_,_]^ + ■ ■ ■ + (a^(i) — l)epy+SQ+...+p._-^+s^_-^+p._i_s.](9.6.12) 
for i = (1, . . . /), and 

(7 — {U-l ^pa+so+-+Pi+si+l -t- ■ ■ • -t- Up(i)fc'po+So+---+p;+S;+p(0 "T 
+ '^1 %0+S0 + -+pi+Si+pW + l + ■ ■ ■ + '^s(i+l)6po+SO + --+Pi+Si+Pi + J + 

+ {{bi^ — l)epy+so+...+p,+s,+p,_|_j+i + • • ■ + {b^Ji^ — l)epjj^_^|j_^..._,_p^_,_^^_^p^^^^^(i) + 
+ {b[ ^ — l)ep|j_,_3jj+..._,_pj^_^^^_p^^^_,^3(i)+i + ■ h (&^(i+i'*) — l)epQ4.so+...+p,+s;+p,_|_j+s;^J (9.6.13) 

Vector V is: 

t; = t;l + ... + t;'+l (9.6.14) 
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Now we can write down Gelfand, Graev and Retakh hypergeometric series corresponding to 
the lattice B and vector v. 



Wj+7iibjH \-ni+ib'+^ 

n.,..i:,.ezj=\ + riib] + ■■■ + m^.bf' + 1) 
Let us compare this series with tau function ( p.6.3| ): 

Fsiv, z) = ci(a, b)gi{z) ■ ■ ■ Q+i(a, b)gi+i{z)T{M, x,yi,..., yi+i) (9.6.16) 

where 

cT\a, h) = r(l - at'^) ■ ■ -1(1 - a^,7.\) )r(l - b?) ■ ■ -1(1 - b%) x 

xr(6!^-^)) ■ ..T{b%i)T{a?) ■ ■■T{a%), z = l,...,l (9.6.17) 



X 



cr+\(a,&) = r(l - a^;^) • ■ -1(1 - a^;i,)Til - a^^) ■ • -1(1 - a^'+V)) 

xr(6f)) ■ ■ ■^(62))^(6^'^) ■ ■■Tib%i) (9.6.18) 

^Po+soH — hpi-i+si-i+pi+i ■ ■ ■ ^pi+siH — hPi-i+si-i+Pi+Si _-|^ i — 2 Z (9 6 19) 



( ^po+soH hpi-i+si-i+l) ■ ■ ■ ( ^po+soH hpi-i+Sj-i+Pi 



(-2:1) ■ ■ ■ {-Zp^} 

yi = ZN-i^i+i, 2 = 2,...,/ + ! (9.6.21) 



^Pl+SlH |-P!+S!+Pi+1 + 1 ■ ■ ■ ^pi+SlH |-Pi+S;+P! + 1+Si+1 



^ ^pi+siH hPi+Si + l) ■ ■ ■ ( 



X (9.6.22) 



Pi+siH hp;+Si+p!+iy 



fi'll^J ^P0+S0 + -+Pi-1+Si-1+1 ^po+3o + ...+p._-^+s-_;^+p{i-l) ^ 

•••9' V 

Po+sa+-+Pi_i+Si_i+p(^-i)+l ^??o+-5o+-+Pz-i+5z-i+??i ^ 

X^ •••7 X 

^^P0+S0 + ---+Pi-1+Si-1+Pi + 1 ^po+sg-l |_p-_-^ + cj-_-^+p.+s(i-l) ^ 

xz^ ' ^ ■■■z^^''^ ^ r9 6 2S^ 

P0+.sn4 kri.- -,-l-Ti.-+.'!(»-l)+1 PO+SoH 



+S0H hPi-l+Si-l+Pi+s('~^'+l PO+^oH hPi-l+Si-l+Pi+Si-1 



for z = (1, . . . Z), and 



pi+si+---+P!+s;+p(')+l ^Pi+si+---+P!+Si+Pi+i ^ 



^^Pl+Sl+---+Pi+Si+Pi + 1+1 ^p^+S;^-) (-p;+S,+p;^;^+s(0 ^ 



X7^ / . . . 7^ =^ ^ ^ fq 941 

^ pi+si+-+p,+s,+p,+i+s(0+l ^Pi+si+-+Pi+Si+Pi+i+Si+i-l l^».u.z,'±; 
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